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ABSTRACT 


A  theory  is  developed  for  predicting  the  time-dependent  sise  and 
shape  of  cracks  in  linearly  viscoelastic,  isotropic  media,  and  its 
validity  is  demonstrated  by  applying  the  theory  to  crack  growth  and 

t 

failure  of  unfilled  and  particulate-filled  polymers.  Starting  with  a 
bounded  solution  for  the  stress  distribution  near  a  crack  tip  in  au 
elastic  body  and  the  extended  correspondence  principle  for  viscoelastic 
media  with  moving  boundaries,  a  simple  equation  is  derived  for 
predicting  instantaneous  crack  tip  velocity  in  terms  of  the 
opening-mode  stress  intensity  factor;  although  the  undamaged  portion 
of  the  continuum  is  assumed  linear,  no  significant  restrictions  are 
placed  on  the  nature  of  the  disintegrating  material  near  the  crack 
tip  and,  therefore,  this  material  may  be  highly  nonlinear, 
rate-dependent,  and  even  discontinuous.  A  further  analysis  is  made  to 
predict  the  time  at  which  a  crack  starts  to  grow,  and  then  some 
explicit  solutions  are  given  for  this  so-called  fracture  initiation 
time,  the  time-dependent  crack  growth,  and  the  time  at  which  gross 
failure  occurs  under  time-varying  applied  forces  and  environmental 
parameters.  Following  a  derivation  of  the  linear  cumulative  damage 
rule,  an  examination  of  its  theoretical  range  of  validity,  and  a 
discussion  of  the  experimental  determination  of  fracture  properties, 
the  theory  is  applied  to  monolithic  and  composite  materials  under 
constant  and  varying  loads.  Some  concluding  remarks  deal  with 
extensions  of  the  theory  to  include  finite  strain  effects,  crack 
growth  in  the  two  shearing  modes  and  in  combined  opening  and  shearing 
modes,  and  adhesive  fracture. 
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1.  Introduction 


/ 


Prediction  of  the  mechanical  response  and  service  life  cf 
viscoelastic  structures  is  very  much  dependent  on  one's  ability  to 
predict  the  time-dependent  growth  of  flaws  or  cracks.  This  growth 
affects  not  only  the  time  of  gross  failure,  but,  in  the  case  of 
composite  materials,  may  also  have  an  appreciable  influence  on  the 
effective  mechanical  behavior  long  before  structural  failure 
occurs  [1]. 

In  this  paper  we  focus  our  attention  on  the  me  ement  of  a  single 
crack  tip  or  front,  and  establish  and  apply  relatively  simple 
equations  for  predicting  the  time  at  which  motion  initiates  (the 
"fracture  initiation  time")  as  well  as  the  instantaneous  tip 
velocity,  both  in  terms  of  the  opening-mode  stress  intensity  factor. 

Interactions  between  two  or  more  cracks  and  the  influence  of 

_  \ 

boundaries  and  obstacles  are  implicitly  taken  into  account  through 
this  stress  intensity  factor.  For  reasons  of  mathematical  simplicity, 
the  material  in  a  small  neighborhood  surrounding  the  crack  tip  is 
divided  into  two  regions:  (i)  a  failure  zone  where  disintegration 
and  eventual  failure  occur,  and  (ii)  a  linearly  viscoelastic, 
macroscopically  homogeneous  and  isotropic  inertialess  continuum. 

This  idealization  and  some  later  mathematical  approximations  are 
motivated  by  certain  characteristics  of  viscoelastic  behavior  which 
are  common  to  many  different  polymeric  materials.  However,  this 
motivation  does  not  mean  the  theory  is  limited  to  polymers;  indeed, 
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it  contains  as  a  special  case  the  classical  theory  of  fracture  for 
elastic  and  elastic-perfectlv  plastic  media  when  the  plastic  zone 
size  is  small  relative  to  crack  length. 

We  shall  not  attempt  to  review  or  even  list  the  many  experimental 
and  theoretical  papers  which  now  exist  on  the  subject  of  crack 
growth  in  viscoelastic  materials.  Instead,  the  reader  is  referred 
to  the  excellent  review  article  by  Knauss  [2],  which  covers  molecular 
processes,  stress  analysis,  crack  propagation  concepts,  heat  evolution 
effects  (which  we  neglect  in  this  paper),  and  environmental  influences. 
Primarily  only  those  studies  which  are  closely  related  to  the  present 
one  are  discussed  here. 

There  already  exist  some  theories  for  predicting  fracture 
initiation  time  and  crack  growth  in  the  opening  mode  in  linear 
viscoelastic  media.  Williams  [3]  used  the  idealized  geometry  of  a 
disintegrating  hollow  sphere  to  establish  the  essential  dependence 
of  fracture  initiation  time  and  initial  crack  velocity  on  external 
loading  history  .  Wnuk  and  Knauss  [4]  deri/ed  an  equation  for 
fracture  initiation  time  of  a  penny-shaped  crack  in  a 
viscoelastic-perf ectly  plastic  solid,  and  gave  solutions  for  a 
Maxwell  material  and  general  linear  viscoelastic  solid  having  a 
time-dependent  and  time- independent  yield  stress,  respectively. 

In  order  to  predict  crack  growth  following  initiation,  Knauss  and 
co-workers  [5-7]  modeled  the  disintegrating  material  near  the  crack 
tip  by  assuming  bi-singular  stress  and  displacement  distributions 
and  using  a  rather  idealized  mechanism  for  calculating  the  transfer 
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of  energy  from  the  linear  viscoelastic  continuum  to  the  intermolecular 
or  interatomic  bonds.  Kostrov  and  Nikitin  [8]  derived  an  equation 
governing  crack  growth  in  a  linearly  viscoelastic-perfectly  plaatic 
solid,  studied  its  consequences  for  the  Griffith  proDlem  of  a  crack 
in  a  large  sheet  under  constant  stress,  and  gave  an  explicit 
solution  for  a  Maxwell  material. 

One  distinguishing  feature  of  existing  theories  on  crack  growth 
is  that  the  failing  material  near  the  crack  tip  is  replaced  by  a 
very  idealized  model.  Another  feature  is  that  solutions  for  the 
time-dependent  crack  size  are  given  numerically  or,  if  analytical, 
are  for  severely  idealized  media  (e.g.,  a  Maxwell  material). 

The  present  study  was  undertaken  primarily  to  remove  these  two 
limitations,  and  is  accomplished  by  introducing  one  new 
assumption:  the  second  derivative  of  the  logarithm  of 
creep  compliance  with  respect  to  logarithm  of  time  is  small  for 
the  linear  viscoelastic  continuum.  To  the  writer's  knowledge,  all 
polymeric  structural  materials  of  technological  interest  satisfy  this 
assumption  over  all  or  most  of  the  time  range  of  variation  of  the 
creep  compli  ~ce.  Even  if  there  are  intervals  over  which  the 
curvature  is  not  small,  the  theory  is  valid  whenever  the  effective 
cime  parameter  (Eq.  (92a))  is  ou^oide  of  these  intervals. 

The  fact  that  the  failure  process  at  the  crack  tip  is  quite 
arbitrary  in  the  present  theory  implies  the  theory  should  be  valid 
for  many  different  materials,  including  polymers  exhibiting  the 
crazing  phenomenon  [9]  and  filled  polymers.  This  point  is  supported 


in  Section  4  by  successfully  predicting  crack  growth  and  failure 
of  an  unfilled  rubber  (Solithane)  and  two  particulate  composites 
(solid  propellant  and  asphaltic  concrete)  under  constant  and 
varying  loads.  Moreover,  because  of  the  simplicity  and  generality 
of  the  theory,  we  believe  it  will  be  useful  for  a  variety  of 
applications,  such  as  the  prediction  of  effective  viscoelastic 
properties  of  composite  media  having  "microstructural  damage"  [1], 
fatigue  crack  growth,  and  the  statistical  distribution  of  failure 
times  of  monolithic  and  composite  materials  in  terms  of 
fundamental  property  and  flaw  size  distributions  [10], 

We  now  turn  to  a  study  of  the  elastic  and  viscoelastic 
distributions  of  stresses  and  displacement  near  the  crack  tip. 

The  nature  of  these  distributions  is  examined  in  some  detail  in 


order  to  provide  a  basis  for  making  certain  later  approximations. 
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2.  Stress  and  Displacement  Distributions  Near  the  Crack  Tip 

< 

The  problem  considered  in  this  Section  is  that  of  predicting 
stresses  and  displacements  in  the  vicinity  of  a  crack  tip  in  which, 
for  simplicity,  the  relative  displacement  between  the  crack  surfaces 
is  taken  normal  to  them;  i.e.,  the  so-called  "opening  mode"  is  assumed. 

A  fixed  (x,y)  coordinate  system  and  the  idealized  tip  geometry 
are  shown  in  Fig.  1;  y  =  0  will  always  correspond  to  the  plane  of  the 
crack,  while  the  y-axis  is  placed  at  any  convenient  location  and 
embedded  in  the  continuum.  In  the  unstrained  state  the  crack  surfaces 
near  the  tip  are  assumed  to  be  planar  and  to  coincide  with  the  x-z 
plane,  where  z  is  perpendicular  to  the  plane  of  the  page.  The 
actual  crack  tip  is  a  straight  or  curved  line  in  space  whose 
intersection  with  the  x-y  plane  (the  plane  of  the  page)  is  defined 
by  the  point  P  (x  =  a,  y  =  0) ;  it  is  assumed  that  the  tip  is 
essentially  straight  in  the  neighborhood  of  P  and  perpendicular  to 
the  x-y  plane.  We  should  add  "neighborhood  of  the  point  P"  refers 
simply  to  the  material  contained  in  a  spherical  volume  centered  at  P 
and  having  dimensions  on  the  order  of  the  failure  zone  width,  a;  P  is 
taken  to  be  a  generic  location  along  the  crack  tip  and,  therefore,  the 
term  "neighborhood  (or  vicinity)  of  the  crack  tip"  will  be  used  to 
mean  the  same  thing  as  "neighborhood  of  P"  throughout  this  paper. 

The  crack  tip  may  be  quite  general  in  shape  as  far  as  our  analysis 
is  concerned.  For  example,  it  may  be  a  closed  curv,e  or  begin  and  end 
on  bounding  surfaces  of  the  continuum;  of  course,  the  first  case 
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corresponds  to  an  "internal"  crack,  while  the  crack  tip  in  the  second 
case  defines  "through"  and  "surface"  cracks. 

We  assume  the  material  outside  of  the  failure  zone  shown  in 
Fig.  1  is  linearly  viscoelastic.  However,  in  general,  the  material 
in  the  zone  may  be  highly  nonlinear  and  viscoelastic;  further,  it 
need  not  be  a  continuum  in  that  it  could  consist  of  very  fine  strands 
or  any  other  disintegrating  material  which  is  representative  of  the 
complex  structure  near  a  crack  tip,  or  the  failure  zone  could  be 
simply  the  region  in  which  significant  attractive  interatomic  forces 
act  between  atoms  in  adjacent  crack  surfaces.  In  practice,  the 
boundary  between  the  undamaged  continuum  and  the  material  in  the 
failure  zone  will  be  rather  diffuse,  and  there  will  be  a  layer  of 
damaged  material  which  extends  into  the  continuum  along  x  <  a^  as 
well  as  around  the  fail, re  zone;  we  assume  this  layer  is  small  enough 

that  its  effect  on  the  prediction  of  the  displacement  v  in  the 
neighborhood  of  the  tip  is  negligible. 

The  location  x  =  a  defines  the  apparent  crack  tip,  which  is  often 
Si  " 

referred  to  in  the  literature  as  the  crack  tip.  However,  in  contrast, 
throughout  this  paper  crack  tip  will  mean  the  tip  at  x  *»  a. 

The  reaction  of  the  failure  zone  on  the  surrounding  continuum 
is  represented  by  a  tensile  stress  distribution,  a^,  as  shown  in 
Fig.  2,  while  o denotes  the  tensile  stress  acting  within  the  linearly 
viscoelastic  material  along  the  crack's  prolongation.  External 
loads  applied  to  the  continuum  are  assumed  to  act  symmetrically  with 
respect  to  the  crack,  so  that  shearing  stresses  along  x  >  a  are  zero 


and  a  pure  crack  opening  mode  results. 


However,  in  reality,  shearing  stresses  act  along  the  continuum 
within  the  failure  zone.  This  may  be  seen  by  referring  to  Fig.  3, 
which  shows  the  forces  acting  on  the  material  in  this  zone;  the  tip 
has  been  magnified  to  bring  out  the  tip  thickness,  d^,  which  is  taken 
to  be  a  characteristic  dimension  defining  an  interatomic  distance 
before  failure  in  the  case  of  a  monolithic  material,  or  an 
interparticle  distance  in  the  case  of  a  particulate  composite.  As 


will  be  seen  later,  the  stress  ax  acting  along  x  >  a  is  equal  to  a 
and,  therefore,  =  °y^i  (Por  un*t  len8th  of  crack  tip).  Thus, 
a  rough  estimate  of  the  average  shearing  stress,  x^,  is 


F  o  d 

_  x  y  i 

f  “  a  2  a 


Consequently,  as  long  as  a  >>  d^,  it  is  reasonable  to  neglect  the 
effect  of  T^,  relative  to  that  of  o^,  on  the  surrounding  linear 
viscoelastic  material;  this  condition  is  consistent  with  that  needed 
to  be  able  to  apply  continuum  mechanics  in  the  first  place  to 
predicting  the  response  of  (assumed)  homogeneous  material  close  to 
the  failure  zone. 

Although  Figs.l  and  2  show  the  crack  surface  for  x  <  a  as 

a 

being  stress-free,  the  theory  allows  for  forces  to  act  directly  on 
the  crack  surface,  as  long  as  they  are  not  applied  in  the  neighborhood 
of  the  tip.  If  external  forces  close  to  the  tip,  the  theory  will 
be  at  least  approximately  valid  when  they  are  small  relative  to  the 
resultant  failure  force.  (It  is  of  interest  to  point  out  that  an 
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I 


important  geological  problem  in  which  there  is  direct  pressure  on 

crack  faces,  viz.,  fracture  of  oil-bearing  rock  strata,  satisfies  this 

condition;  there  is  always  an  unwetted  portion  near  the  crack  tips  [11].) 

Finally,  it  is  assumed  that  the  neighborhood  of  the  point  P  in 

Fig.  1  is  in  a  state  of  plane  strain.  This  condition  is  (approximately) 

satisfied  if  (i)  the  distance  between  P  and  bounding  surfaces  of 

the  continuum  and  other  cracks  is  large  compared  to  a;  (ii)  the  radius 

of  curvature  of  the  crack  tip  at  P  is  much  larger  than  a;  and 

(iii)  the  value  of  outside  the  neighborhood  of  P  is  small  relative 

to  0^  inside  the  neighborhood.  Of  course,  if  one  is  interested  in 

the  problem  of  "through  crack"  in  a  thin  sheet  in  which  a  is  large 

compared  to  sheet  thickness,  the  plane  strain  results  are  converted  to 

this  plane  stress  situation  by  means  of  a  simple  change  in  the  solution; 

2 

i.e.,  v  in  Eq.  (7)  is  dropped. 


Linear  Elastic  Analysis 


Stress  and  displacement  distributions  along  y  =  0  in  the 

neighborhood  of  the  crack  tip  are  recorded  here  for  later  extension  to 

viscoelasticity.  We  first  consider  stresses  o°  and  o°,  and  displacement 

x  y 

v  (say)  due  to  loads  acting  in  and  on  a  linear  elastic  continuum 
having  a  crack  tip  at  x  *  a  and  with  =  0.  Based  on  the  results 


given  by  Williams  [12],  we  may  write* 


—  H(51)[l  +  0(£  /B)] 
1 


*A  term  of  order  p  is  written  as  0(p)  for  which,  by  definition, 
lim  [0(p)/p]  ■  finite,  non-zero  quantity. 

p  0 
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N 

o°  *  H(L)[1  +  0<£  /{$)] 

y  ^7 


while  the  displacement  of  the  top  crack  face  is  (for  plane  strain) 


v°  =  Ce  Nq  /K  H(0[1  +  0(5/B)] 


where  (see  Fig.  2) 


5^  =  x  -  a  ,  £  =  a  -  x 


Also,  H(p)  is  the  unit-step  function: 


H(p)  = 


0  ,  p  <  0 
1  ,  p  >  0 


4(1  -  vi 


where  E  is  Young's  modulus  and  v  is  Poisson's  ratio.  The  coefficient 
Nq  is  called  the  stress  intensity  factor  by  Barenblatt  [11];  this 
factor  is  simply  related  to  another  commonly  used  symbol,  Kj.,  for 
the  opening-mode  stress  intensity  factor 


Kt  =  /2tT  N 
I  o 


In  general,  Nq  is  a  functior  of  Young's  modulus  and  Poisson's  ratio 
as  well  as  parameters  which  define  the  geometry  of  the  continuum 
(e.g.,  crack  length  and  location  and  overall  dimensions  of  the  body) 
and  is  a  linear  function  of  the  applied  loads  and  displacements. 


The  quantity  $  represents  the  distance  between  the  point  F 
(Fig.  1)  and  the  closest  geometric  feature  (e.g.,  free  surface  or 
opposite  crack  tip)  or  nearest  point  of  load  application,  depending 
on  which  is  smallest.  By  restricting  £  and  £^  to  the  neighborhood 
of  F  and  assuming 

a/$  «  1  (9) 

the  terms  of  order  £,/6  and.  £/B  can  be  neglected  since 

a 

°<r> '  0(t  h  0<f)-0<tr)-  (l0> 


Hence,  close  to  the  crack  tip 


a 


o 

X 


(11) 


V°  -  c  N  /l  H(£)  . 

e  o 


(12) 


Stresses  and  displacement  in  the  linear  solid  due  to  the  failure 

f 

y’ 


stress,  o  ,  acting  alone  are  denoted  by  o^,  o^,  and  v^.  From 
i  x  y 


Barenblatt  [11] ,  for  £^  >  0, 


f  f 

o*o 
x  y 


if 


°  ^  of(£)  d£ 

e  +  5, 


(13) 


and  the  displacement  of  the  top  crack  face  is 
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f 

V 


-■57  »«>«• 


o 


ofa')  ln[- 


/c  +  si 

Sc  -  SI 


]  dC7 


(14a) 


where  &  denotes  "real  part  of";  inasmuch  as 

ln(-l)  «  i  it  (2k  +  1),  k  «  0,1,2,...,  the  last  equation  can  be  written 
in  real  notation  as 


f 

v 


57  ’«» 


■a 


of(5’)  In 


SI7  +  /[ 

Sc  - 


d£' 


(14b) 


On  physical  grounds  it  is  assumed  that  o»  is  everywhere  finite, 
which  assures  uniform  convergence  of  the  above  integrals,  Eqs.  (13)  and 
(14),  for  £  >  0  and  £  >  0.  Although  the  failure  stress  is  not 

restricted  to  being  continuous  for  all  £,  we  do  make  the  physically 
plausible  assumption  that  it  is  continuous  to  the  left  at  £  =  0;  viz. 

lim  0  (O  -  of  (0)  .  (15) 

4  0+  f  f 

As  will  be  seen  later,  existence  of  the  limit  in  Eq.  (15)  together 
with  the  assumption  that  the  stresses  in  the  linear  continuum  are 
finite  establishes  continuity  of  the  tensile  stress  in  the  y-direction 
at  y  =  (  =  0;  even  if  Eq.  (15)  is  not  mec,  it  can  be  shown  that  the 
stress  immediately  to  the  right  of  the  tip,  o^(0+),  still  will  be 
equal  to  the  failure  stress  on  the  left  side  of  the  tip,  o^(0+). 

Of  particular  interest  to  us  is  the  behavior  of  stresses,  Eq.  (13), 
and  displacement,  Eq .  (14),  close  to  the  crack  tip;  viz.,  4^/«  <<  1  and 
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£/a  <<  1.  The  stresses  become  [11] 

ra 

f  f  1  wo 

°  -  o  -  -  d£  +  o  (0)  +  o(d^)  .  (16) 

y  J0  »T  1  1 

Although  not  stated  in  [11],  the  order  of  the  remainder  term  in 
Eq.  (16),  as  well  as  the  remainder  shown  later  in  displacement 
Eq.  (25),  is  correct  only  if  satisfies  the  following  condition: 

lim  { [a  (O  -  af(0)j/£1/2}  -  0  (17) 

which  will  be  established  and  discussed  later  in  this  subsection.  The 
resultant  stresses  in  the  linear  continuum  due  the  combined  action 
of  the  failure  stress  distribution  and  the  externally  and  internally 
applied  loads  are  obtained  by  adding  Eq.  (11)  and  (16).  These 
stresses  will  be  finite  at  the  tip  if  and  only  if 


A  more  complete  representation  is  obtained  by  adding  Eqs.  (11)  and 
(13)  and  then  using  Eq.  (18): 
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Since  £  <  1, 


I 


2 


(24b) 


which  will  prove  useful  in  later  considerations  in  order  to  establish 
a  lower  bound  on  the  maximum  strength,  o^,  from  experimental  data. 

The  normalized  form,  Eq.  (22),  is  a  very  Important  result.  It 
will  enable  us  to  calculate  the  length  of  the  failure  zone,  a,  as  a 
function  of  time  and  loading. 

The  displacement  of  the  top  crack  face  near  the  tip  due  to 
Of  acting  alone,  as  given  by  Barenblatt  [11],  is 


f 


v 


-  ~  H(C) 


; -  +  Av*  H(0 


(25) 


f  3/2 

where  Av  ■  0(5  ).  In  order  to  predict  time-dependent  crack  growth 

in  viscoelastic  materials  we  need  more  explicit  representations  of 

Av^;  for  one  purpose  it  will  be  necessary  to  evaluate  the 
3/2 

coefficient  of  5  and  establish  the  order  of  the  remaining  terms. 
This  will  be  accomplished  by  first  equating  Eqs.  (14b)  and  (25),  which 
yields 


Av1, 


C 

e 

2ir 


•a 

’ 

" 

ofa')  ■ 

2(^1 1/2  - 

In 

Jv  +  Ji 

f 

o 

V 

» 

Jv  -Jl 

- 

(26) 


v 


It  is  important  to  observe  that  the  resultant  displacement, 
v°  +  v^,  due  to  applied  loads  and  to  o^  is  simply  related  to  Av^; 
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viz. , 


v  »  Av  H(£) 


which  follows  from  Eqs.  (12),  (14b)  and  (18). 

Determination  of  the  behavior  of  Avf  from  Eq.  (26)  requires  that 


we  express  a  in  the  form 


°f  (£)  *  of(0)  +  A  a  (O 


In  view  of  Eq.  (15), 


lira  Ao  (5)  =  0  . 

5  -►  0+  f 


Upon  substituting  Eq.  (28)  into  Eq.  (26)  and  evaluating  the  integral 


containing  of(0)  there  results 


*vf .  £  tsw  +  „«*«>,  +  „ 

*a 


where 


*  \  Jo  /,1_\  1/2 


A!  =  Aof  (5 ' )  ^2(|t)1/2  -  in  Ml 


1  -  /?/?’ 


d?'  .  (31) 


The  term  in  curly  brackets  for  0  <  C/5'  «  1  is  approximately 


2 (It)  1/2  -  In  L  .  2  £_  3/2  _  2  £_  5/2 

5  l  -  /^yv  [  3  5  V}  (32) 

j 


Without  any  approximation,  Eq.  (31)  is  now  written  as 


16 


ra 


at  -  e  r3/ 2 
AI  "  ’  37  C 


Aof(£‘) 

(c’)3/2 


d£' 


ca 


Aof(£')  |2(|r)1/2  -  m 


1  +  /c7T 


1  -  JUT 


+  1  (i-)3/2 

3  V; 


d£’ 


(33) 


If  Eq.  (17)  is  satisfied  tiie  first  integral  in  Eq.  (33)  obviously 
converges;  also,  the  second  integral  is  easily  shown  to  be  convergent. 

This  condition  means,  of  course,  that  Ao^  must  approach  zero  faster 

1/2  3/2 

than  £  as  C  0.  If  this  approach  to  zero  is  faster  than  £  ; 

viz. , 

lim  [A of(£)/£3/2]  -  0  (34) 

€  0+  f 

then  it  follows  from  Eqs.  (32)  and  (33)  that  the  second  integral  in 

5/2 

Eq.  (33)  is  of  order  £  or  higher.  In  contrast,  if  Eq.  (17)  is 

satisfied,  but  not  Eq.  (34),  then  the  order  of  the  second  integral  ir. 

3/2  5/2 

Eq.  (33)  is  found  to  be  between  £  and  £  .  Moreover,  it  is  shown 

in  Appendix  A  that  if  Eq.  (17)  is  not  met,  the  order  of  AI,  Eq.  (31), 
3/2 

is  between  £  and  £ 

Considering  the  complexity  of  the  failure  process  at  the  crack 
tip,  and  the  absence  of  information  on  the  variation  of  o^(£)  for 
most  materials  [2],  it  will  be  assumed  throughout  the  remainder  of 
this  paper  that  Eq.  (17)  holds,  unless  stated  otherwise.  Note  that 


if  is  constant  (eg.,  when  the  failure  zone  is  actually  the  zone 
of  yielding  in  a  metal  having  zero  strain  hardening)  not  only  will 
Eq.  (34)  be  satisfied  but  additionally  41  *  0. 

Let  us  now  integrate  the  first  integral  in  Eq.  (33)  by  parts, 
use  definition  Eq.  (28)  for  Aa^,  replace  the  second  integral  by  an 
order  statement,  and  then  substitute  the  result  into  Eq.  (30);  thus, 

Avf  .  _£  ^.U3/2  +  0(£3/2  +  p)  (35: 


where  0  <  p  <  1,  and  0£(a“)  is  the  failure  stress  for  £  just  to  the 
right  of  a;  of  course,  if  is  continuous  at  £  *  a  there  is  no  need 
to  distinguish  between  £  -  a  and  £  *  a.  Equation  (35)  can  be  written 
-lore  compactly  in  the  form 

/vf  =  3/2  ,3/2  +  p. 

Av  TT  dC  7=  e  +  °a  }  *  (36) 

■  o 

It  is  easily  verified  that  Eq.  (36)  is  equivalent  to  Eq.  (35);  when 
is  discontinuous  at  £  *  a  we  may  interpret  do^/d£'  by  means  of  a 
Dirac  delta  function  or,  more  directly,  by  using  the  relation 


Z1  dil 
dC  /£T 


cT  .  or  , 

i hiV..  tl 

JFT  +  <U 

o  ^  }a~ 


_i  d£l 

d5'  /c 


f  d£'  1  sf. 


oAoT)]  .  (37) 


i  lii  rmr— v~‘-1 


The  total  vertical  displacement  of  the  upper  crack  face  due  to 
external  and  internal  loading  and  to  the  failure  stresses  is  given  by 


Eq.  (27).  After  writing  the  integral  in  Eq.  (36)  in  terms  of  the 
normalized  variables,  Eqs.  (20)  and  (21),  there  results 


2C  a 

3—  -p  I253/2H(S)  +  0(C3/2  +  P)  H(0 


(38) 


where  0  <  p  <  1  and 
1 

df  dp 

dn  r~ 
/  0 


(39) 


The  upper  limit  in  Eq.  (39)  is  to  be  interpreted  as  1+  when  f  is 
discontinuous  at  n  1  1;  for  example,  if  (and  therefore  f)  is 
constant  over  0  <  n  <  1, 


I2  -  -1  •  (40) 

Also,  p  *  1  for  this  case  since  AI  =  0. 

3/2 

The  fact  that  v  =  0 ( C  )  according  to  Eq.  (38),  the  cross-section 

of  the  crack  near  its  tip  (see  Fig.  1)  is  cusp-shaped.  This  observation 
was  made  by  Barenblatt  [11].  However,  as  we  have  shown,  this  shape  is 
predicted  because  Eq.  (17)  is  assumed  to  be  satisfied.  On  the  other 
hand  when  this  limit  condition  if  not  met,  a  study  of  AI,  Eq.  (31), 
in  Appendix  A  reveals  that  the  cros  section  could  approach  a 
straight-sided  wedge. 

Finally,  we  should  point  out  that  the  order  of  the  stresses  near 
the  tip  can  be  established  in  a  manner  that  is  entirely  analogous  to 
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the  above  one  for  displacement  by  starting  with  Eq.  (19b);  it  is  found, 

for  example,  that  the  remainder  term  in  Eq.  (19a)  is  indeed  of  order 

1/2 

5^  (or  higher  order  when  I^,  Eq.  (39),  vanishes)  if  Eq.  (17)  is 
obeyed.  That  the  order  of  this  term  could  range  between  £°  and  $ 
when  this  limit  condition  does  not  hold  is  shown  in  Appendix  A. 

The  order  of  the  stresses,  however,  does  not  have  a  direct  bearing  on 
the  subsequent  viscoelastic  analysis. 

Linear  Viscoelastic  Analysis 

With  stationary  cracks  the  elastic  stress  and  displacement 
distributions  given  in  the  previous  subsection  can  be  easily  generalized 
to  viscoelastic  solutions  by  means  of  the  classical  correspondence 
principle  plus  Laplace  transform  inversion  [13]. 

Under  somewhat  restrictive  conditions,  the  extended  correspondence 
principle  established  by  Graham  [14]  can  be  used  with  moving  cracks. 

The  specific  restrictions  as  they  apply  to  the  problem  at  hand  are 
(a)  the  crack  cannot  decrease  in  size  (da/dt  >  0);  (b)  the  elastic 
stress  normal  to  the  surface  of  crack  prolongation  (o^)  must  be 
independent  of  E  and  v;  and  (c)  any  dependence  of  the  elastic 
displacement  v  along  the  crack  face  on  E  and  v  must  be  in  the 
form  of  a  separate  factor  (i.e.,  v  =  fntn(E,v)  •  fntn(£))*  According 
to  the  theory  in  [14],  the  above  restrictions  apply  to  the  entire 
crack  and  its  extension  up  to  the  maximum  time  of  interest,  and  not^ 
just  the  neighborhood  of  the  instantaneous  tip. 

With  reference  to  restriction  (b) ,  we  note  that  as  long  as  all 
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boundary  conditions  are  on  stress,  dimensional  analysis  implies  the 
internal  stresses  will  be  independent  of  Young's  modulus,  E;  they 
may,  in  general,  depend  on  Poisson's  ratio,  v.  Even  this  dependence 
on  v  disappears  in  many  two-  and  three-dimensional  problems. 

Indeed,  there  are  two  important  classes  of  plane  stress  and  plane 
strain  problems  for  which  the  in-plane  stresses  will  always  be 
independent  of  v  [15]:  (i)  simply-connected  bodies  and 

(ii)  multiply-connected  bodies  for  which  the  resultant  force  on  each 
boundary  vanishes.  It  should  be  added  that  if  the  viscoelastic 
material  of  interest  has  a  constant  Poisson's  ratio  (which  is  at 
least  approximately  true  in  many  situations  involving  polymeric 
materials  [  1  ]),  then  dependence  of  the  elastic  solution  on  Poisson's 
ratio  does  not  invalidate  the  extended  correspondence  principle; 
this  last  point  is  easily  verified  by  examining  Graham's  theory  [14]. 
As  one  further  observation,  we  note  that  if  some  or  all  boundary 
conditions  are  on  displacement,  the  problem  often  can  be  recast  in 
terms  of  applied  forces;  one  simply  calculates  the  resultant  forces 
needed  to  maintain  the  specified  displacements,  and  then  views  these 
forces  as  being  specified  quantities.  It  will  be  assumed  that  such  a 
recasting  has  been  accomplished,  when  appropriate,  in  all  of  the 
following  analyses;  see,  e.g.,  Eqs.  (128)  and  (129). 

Under  the  conditions  for  which  the  extended  correspondence 
principle  is  applicable,  the  stresses  in  a  linearly  viscoelastic  body 
are  the  same  as  those  given  in  the  preceding  subsection  for  elastic 
media.  The  only  stress  result  which  will  be  needed  for  predicting 
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crack  growth  Is  Eq.  (22)  relating  the  stress  Intensity  factor,  Nq, 
tc  the  failure  stress  and  size  of  the  failure  zone.  Recall  that, 
by  definition,  Nq  is  the  stress  intensity  factor  for  a  crack  with 
of  =  0.  This  factor  depends,  at  most,  on  parameters  defining  the 
geometry  of  the  body  (including  instantaneous  crack  size  and  location) 
and  on  the  applied  stresses;  in  view  of  the  earlier  comment,  Poisson's 
ratio  dependence  is  not  precluded  as  long  as  it  is  constant. 

Turning  to  the  prediction  of  displacement  v  along  y  ■  0, 
we  observe  that  this  displacement  close  to  the  crack  tip,  Eqs.  (26) 
and  (27),  meets  restriction  (c) .  For  a  stationary  or  growing  crack, 
the  extended  correspondence  principle  yields  for  t  >  t.: 


'•h  Cy(t  ■  ’*  It  '  "£«•>  [2(K/2 


in  ]  d£'  ■  dx 

Jv  - 


The  variable  of  integration  t  is  a  generic  value  of  time  which  ranges 
from  the  time,  t^,  when  the  crack  tip  first  reaches  a  point  x,  to  the 
current  time  t;  this  time  integration  is  to  be  performed  at  constant  x, 
and  therefore  for  moving  cracks  5  must  first  be  expressed  in  terms  of 
x  and  t;  viz.,  from  Eq.  (5): 


£  =  £(jc,t)  -  a(x)  -  x 


In  view  of  the  above  definition  of  t^, 
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^(x.t^)  -  a(t^)  -  x  -  0  (43) 

and  since  decreasing  crack  size  is  not  allowed, 

£(x,t)  >  0  ;  t  >  t1  .  (44) 


Also,  although  not  explicitly  shown,  the  length  of  the  failure 
zone  a  and  the  failure  stress  may  depend  on  time,  x. 

If  the  crack  is  stationary  we  may  set  t^  ■  0,  corresponding  to 
the  time  at  which  loads  are  first  applied  to  the  body;  Eq.  (41)  then 
applies  for  x  <  a. 

The  function  Cy(t)  is  defined  by  the  equation 

s  C  i  W  •  _fr>2)  (45 

V  sE 

in  which  the  bar  denotes  the  Laplace  transform  (LT)  of  a  function  of 
time;  e.g., 


e  St  q(t)  dt 


(46) 


‘  o 

Also,  E  is  the  LT  of  the  uniaxial  relaxation  modulus  and  v  is  the  LT 
of  the  Poisson's  ratio  for  a  uniaxial  relaxation  test.  Equation  (45) 
can  be  easily  established  by  comparing  elastic  stress-strain  equations 
with  the  LT  of  those  for  a  viscoelastic  material  [13].  For  some 
problems,  it  may  be  desirable  to  express  C^(t)  in  terms  of  the  uniaxial 


creep  compliance  D(t),  where 


sD  * 


1_ 

sE 


(47) 


If  Poisson's  ratio  Is  constant  with  the  value  v,  then  sv  ■  v  and 
we  can  combine  Eqs.  (45)  and  (47)  to  find,  by  LT  Inversion, 

Cy(t)  -  4(1  -  v2)  D(t)  .  (48) 

Even  when  Poisson’s  ratio  is  not  constant,  Eq.  (48)  with  v  ■  v(t)  is 
expected  to  be  a  good  approximation  [  1  ];  this  result,  which  is  based 
on  the  quasi-elastic  method,  can  be  expressed  just  as  well  using  the 
Poisson's  ratio  measured  in  a  uniaxial  creep  test,  vc(t),  since 
typically  vc(t)  =  v(t)  (  1  ] . 

Recall  that  Eq.  (41)  is  based  on  an  elastic  solution  which  is 

valid  only  in  the  neighborhood  of  the  crack  tip.  Therefore  one 

cannot  expect  this  viscoelastic  solution  to  be  valid  for  values  of 

£  much  greater  than  a;  this  limitation  does  net  restrict  the  subsequent 

fracture  analysis  as  we  will  need  the  displacement  only  in  the  interval 

0  <  £  <  a.  An  additional  result  that  will  be  needed  is  the  viscoelastic 
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version  of  Eq.  (38)  in  which  only  the  term  of  order  £  is  retained: 


2_ 

3ir 


Cv(t 


,  3 
T)  17 


£  **  53/2 


dt 


(49) 


where  t  >  t^.  All  of  the  comments  following  Eq.  (41)  apply  as  well 
to  Eq.  (49)  except,  strictly  speaking,  this  latter  equation  is 
restricted  to  the  range  0  <  £  <<  a;  in  Appendix  C  we  demonstrate  that, 
in  reality,  Eq.  (49)  represents  quite  a  good  approximation  for  the 
purpose  of  calculating  work  input  to  the  failure  zone. 

Finally,  it  should  be  mentioned  that  restriction  (c)  of  Graham’s 
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extended  correspondence  principle,  which  was  stated  early  in  this 
subsection,  is  overly  restrictive.  Namely,  the  product  dependence  of 
the  displacement  on  material  properties  is  not  essential,  as  we  shall 
presently  show.  This  latter  observation  is  particularly  important 
to  us  because  the  elastic  displacement  in  Eqs.  (26)  and  (27),  which 
has  the  product  dependence,  applies  only  to  the  neighborhood  of  the 
crack  tip;  whether  or  not  this  dependence  exists  far  away  from  the 
tip  probably  depends  on  the  specific  problem  considered,  and  we  do 
not  want  to  restrict  the  viscoelastic  displacement  solution 
unnecessarily.  Now,  following  the  reasoning  in  Graham's .paper  [14]  and 

f 

assuming  restrictions  (a)  and  (b)  are  met,  we  conclude  that  the  LT  of 
the  viscoelastic  displacement  anywhere  along  the  entire  crack  plane  is 
v  *  v  [sE,  sv,  x,  s]  where  v  is  the  LT  of  the  elastic  displacement 
in  which  Young's  modulus  and  Poisson's  ratio  have  been  replaced  by 
sE  and  sv,  respectively.  Explicitly, 


—  —  —  g  £ 

v^lsE,  sv,  x,  t]  e  dt 


(50) 


where  we  have  used  the  fact  that  the  elastic  displacement  at  the 
point  x  along  the  line  of  crack  prolongation  is  zero  until  time  t^. 
Change  the  variable  of  integration  to  p  =  t  -  t  : 


r8ti 


v^tsE,  sv,  x. 


P  +  t1] 


-SP 


dp 


(51) 


Denoting  the  inverse  by },  Eq.  (51)  yields 


v=^-M 


-sp  , 

v  e  dp 
e 


H(t  -  tx) 


(52) 


which  shows  Chat  this  displacement  correctly  vanishes  for  t  <  t 
While  there  is  no  need,  in  principle,  to  assume  a  special  form  for 
the  dependence  of  ve  on  sE  and  so,  the  Inversion  is  greatly  facilitated 
if  this  displacement  consists  of  a  sum  of  terms,  or  a  single  term,  in 
which  the  dependence  on  E  and  v  factors.  Specifically,  for  th. 
problem  at  hand  we  can  write 

f“  — 

-sp  8CV  —  _ 

»e  e  dp  -  —  „  +  4ve  (53) 

o  ' 

where  q  is  the  LT  of  the  integral  in  Eq.  (26),  in  which  ?  =  a(p  +  - 

also,  Avfi  is  the  LT  of  the  higher  order  terms  in  Eq.  (4). 

Substituting  Eq.  (53)  into  (52),  neglecting  Avfi,  using  the  convolution 
rule  for  the  inversion  of  a  product  of  transforms  and  replacing  p  by 
t  -  t1  yields  Eq.  (41) . 

3.  Analysis  of  Crack  Motion 

Equations  for  predicting  time-dependent  crack  extension  will  be 
developed  in  this  section;  the  preceding  analysis  plus  an  energy 
criterion  for  failure  provide  the  basis  for  these  predictions. 

However,  we  believe  it  will  be  helpful  to  first  discuss  briefly  the 
nistory  of  events  which  occurs  in  the  failure  zone  from  the  standpoint 
of  the  simplified  model  illustrated  in  Fig.  1. 

Suppose  an  Initially  unstressed  md  unstrained  body  Is  subjected 
to  gradually  Increasing  load(s).  The  stress  intensity  factor  in 
Ed-  (22)  will  then  increase  In  time;  at  any  instant  the  length  of  the 


26 


failure  zone  must  satisfy  this  equation,  even  before  the  left  end 

(S  *  a)  of  the  zone  starts  to  move.  (If,  for  example,  the  maximum 

failure  stress,  o  ,  and  the  "shape"  of  the  failure  stress  distribution, 
m 

as  defined  by  f  in  Eq.  (21),  remain  constant,  we  see  a,  and 

2 

the  efore  a,  will  grow  in  proportion  to  Nq.)  Eventually,  with 
increasing  load,  a  smai L  column  of  material  at  the  left  end  of  the 
failure  zone  will  rupture  and  the  entire  failure  zone  will  start  to 
move  to  the  right;  the  time  at  which  this  movement  initiates  will  be 
called  "fracture  initiation  time,"  and  denoted  by  t^.  In  many  cases 
the  time  t^  will  be  extremely  small  compared  to  the  total  time  for 
complete  fracture  of  a  body  [2]  and,  therefore,  may  not  be  of  direct 
practical  use;  prediction  of  the  eventual  crack  growth  behavior  would 
then  be  of  primary  interest  for  establishing  gross  fracture  times. 

In  developing  the  equations  governing  crack  growth  under  various 
loading  conditions,  it  is  desirable  to  divide  the  analysis  into  two 
parts:  (i)  continuous  crack  growth  corresponding  to  times  t  >  t^ 
and  (ii)  initial  and  intermittent  crack  growth  corresponding  to  times 
0  <  t  <  t^  and  to  periods  during  rfhich  the  applied  loading  drops  below 
that  necessary  for  continuous  crack  growth.  It  should  be  noted  that  the 
loading  may  vary  in  time  in  both  parts  (t)  and  (ii) ;  however,  it  is 
assumed  in  the  first  part  that  the  loading  is  at  least  great  enough 
to  sustain  continuous  motion  of  the  entire  failure  zone  in  the  direction 
of  +x . 


The  same  fracture  criterion  will  be  used  in  both  parts. 
Specifically,  a  column  of  material  with  cross-sectional  area  dA 
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in  the  failure  zone  0  <  y  <  v  is  said  to  fail  when  the  total  work  done 
on  it  by  the  adjacent  linear  continuum  reaches  the  value  T  dA. 

(There  will  be  no  need  to  separately  identify  the  failure  of  perfectly 
brittle  materials.  In  this  case  T  dA  could  be  interpreted  as  the 
mechanical  work  necessary  to  completely  separate  adjacent  surfaces 
with  area  dA,  and  T  would  be  the  so-called  "cohesive  surface  energy"; 
as  indicated  previously,  the  failure  zone  for  this  case  is  simply  the 
zone  in  which  significant  forces  act  between  atoms  in  adjacent  crack 
surfaces.) 

In  general,  the  quantity  T  depends  on  the  entire  history  of 
the  force,  a ^  dA,  applied  to  a  given  column  of  material.  Fortunately, 
however,  there  are  many  viscoelastic  materials  for  which  experimental 
and  theoretical  considerations  indicate  that  T  is  at  least 
approximately  constant  [2].  The  quantity  T,  whether  history-dependent 
or  not,  will  be  called  fracture  energy  in  the  remainder  of  this  paper. 
Equations  which  relate  7.’  to  motion  of  the  surrounding  continuum  are 
given  in  the  next  s_bsection.  We  then  take  up  parts  (i)  and  (ii)  in 
separate  subsections. 

Calculation  of  Work  Done  on  Failure  Zone 

The  energy  T  dA  is  equal  to  the  work  done  by  the  upper  portion 

(y  >  0)  of  the  continuum  against  the  force,  dA,  at  a  fixed  generic 

location,  x,  as  it  suffers  the  total  vertical  displacement,  v^;  this 

displacement  is  defined  by  the  condition  that  ac  -  0  when  v  >  v  .  Thus, 

i  ra 

the  work  input  per  unit  of  new  area  is 


I 

i 


s 
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v 

r  m 


0^  dv 


(54) 


or,  equivalently, 


h 


8v  , 

°f  3t  dt 


(55) 


where  vertical  motion  begins  at  t  and  reaches  v  at  t.;  note  that  the 

J-  m2 

time  derivative  is  taken  with  x  held  constant. 

Another,  more  useful  form  with  continuous  crack  propagation  is 
obtained  by  considering  the  displacement  in  Eq.  (54)  to  be  a  function 
of  the  independent  variables  (x,C),  instead  of  (x,t),  where  by 
definition 


C  *  £(x,t)  =  a(t)  -  x  .  (56) 

It  is  assumed  here  that  ^  <  t  *  t2  and  (a  -  a)  <  x  <  a.  We  select 
x  such  that  it  is  at  the  crack  tip  when  t  -  t^  note  that  ^(x,^)  -  0, 
and(a(t2)  -  x)  is  the  length  of  the  zone  of  failure,  a. 

Equation  (54)  now  becomes 


r 


r  a 


3v 
Jf  3C 


d£ 


(57) 


in  which  v  *  v(x,£)  and  da/dt  >  0. 

The  above  change  of  independent  variables,  (x,t)  (x,£),  may 

be  viewed  as  a  coordinate  transformation;  viz., 


- - - ...  - 


—  •» — 1- - 
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x'  *  x  (58a) 

K  -  a(t)  -  x  .  (58b) 

The  Jacobian  determinant  is  found  to  be  da/dt,  and  therefore  Eq.  (58) 
defines  an  admissible  (one-to-one)  transformation  [13]  as  long  as  the 
crack  velocity  does  not  vanish. 

Prediction  of  Continuous  Crack  Growth 

The  governing  equations  for  crack  growth  will  be  derived  by 
bringing  together  the  finite  stress  condition,  Eq.  (22),  the 
viscoelastic  displacement,  Eq.  (41)  (and  Eq.  (49)),  and  the  fracture 
energy,  Eq.  (57).  Even  if  a  constitutive  equation  for  the  failure  zone 
were  known,  the  task  of  finding  the  solution  a  =  a(t)  to  this  set  of 
nonlinear  integro-dif ferential  equations  would  be  very  formidable,  to 
say  the  least.  There  are,  however,  two  cases  when  the  equations 
greatly  simplify.  / 

The  first  one  is  for  a  material  with  infinite  strength  (o  -+■  °°)  . 

/  m 

It  is  shown  in  Appendix  B  that  the  exact  result  is  the  classical 
condition  for  critical  stress  intensity  factor  of  a  brittle  material  in 
plane  strain  : 
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f  r  E<Q) 

i»[l  -  v2(0)J 


where  E(0)  and  v(0)  are  the  initial  values  of  relaxation  modulus  and 
Poisson's  ratio.  The  stress  predicted  by  this  equation  is  unrealistic- 
ally  high  compared  to  the  values  at  which  cracks  can  be  made  to 
propagate  in  viscoelastic  materials  [  2  },  and  moreover  the  crack 
velocity  is  undefined.  Except  for  representing  a  limiting  situation, 
this  case  of  infinite  strength  at  the  tip,  therefore,  is  of  little 
practical  interest. 

The  second  case  is  based  on  a  certain  property  of  the  creep 
compliance  which  exists  for  most  polymeric  materials.  Specifically, 
a  double-logarithmic  plot  of  creep  compliance  (e.g.,  in  simple  shear 
or  uniaxial  tension)  has  small  curvature  over  most,  if  not  all,  of 
its  range  of  variation  [e.g.,  1  ,  16  ] ;  the  specific  meaning  of 
"small"  will  be  brought  out  later  in  the  analysis.  The  creep 
compliance,  Cv(t),  which  is  defined  by  Eq.  (45),  can  be  expected  to 
have  small  curvature  when  D(t)  has  this  property  since  Poisson's 
ratio  typically  is  a  very  weak  function  of  time;  see  also  Eq.  (48). 

As  one  example,  the  creep  compliance  in  the  glass-to-rubber  transition 
range  for  Solithane  50/50  [6]  under  uniaxial  tension  is  shown  in 
Fig.  4;  this  material  is  an  unfilled,  crosslinked,  amorphous  polymer. 
The  compliance  of  particulate  f illed-polymers  (e.g.,  solid  propellant) 
in  the  glass-to-rubber  transition  range  and  the  compliance  of  plastics 


typically  vary  over  a  greater  range  of  log  t  and  have  even  less 
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curvature  than  seen  in  Fig.  4.  In  view  of  the  generality  of  this 

behavior  and  the  analytical  simplification  it  produces,  we  shall 

concern  ourselves  with  this  case  in  the  remainder  of  this  subsection. 

As  noted  previously,  a  and  may  be  functions  of  time,  although 

this  dependence  is  not  explicitly  shown  in  Eqs.  (41)  and  (49). 

However,  it  will  be  assumed  that  the  time  it  takes  for  the  crack  to 

propagate  an  amount  a  is  small  enough  that  a  and  are  essentially 

timewise  constant  during  this  brief  period,  denoted  as  t^  -  t. . 

Note  that  Eq.  (18)  (or  Eq.  (22))  implies  N  must  be  essentially 

o 

constant  during  this  same  period;  nevertheless,  it  is  important  to 
realize  that  may  vary  either  weakly  or  strongly  with  time  during 
the  total  period  of  crack  growth.  Additionally,  we  assume  a(t)  can 
be  linearized  during  the  short  period  t^  -  t^;  thus  with  0  <  £  <  a, 

a(t)  ~  a(t^)  +  (t  -  t^)a  (60a) 

and  using  Eq.  (42), 

£(x,t)  =  a(t)  -  x  =  (t  -  tx)a  (60b) 

£(x,t)  =  a(t)  -  x  =  (t  -  t^a  (60c) 

where  a(t^)  =  x  and  a  is  the  tip  velocity  at  the  generic  time  t^. 

We  consider  t^  to  be  a  function  of  x  since  it  is  the  time  at  which 


the  crack  tip  reaches  the  fixed  location  x. 
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Pv  =  p/t 
L  =  log  t 


l  -  log (1  -  pv) 


Cv(L)  i  Cv(t) 


(64) 


Note  that  Cv(L)  is  simply  the  creep  compliance  expressed  in  terms  of 
logarithmic  time.  Equation  (63)  becomes 


fo 


C(  ,(t)  =  j  (In  10) 


C  (L  +  l]  -'10l  [1  -  105']1/2\  tl. 
v  l  1  / 


(65) 


The  quantity  in  curly  brackets  will  be  called  a  weighting  function 
and  denoted  w^^I  plotted  in  Fig.  5  over  its  major  range  of 

variation.  Recognizing  that  tne  creep  compliance  is  an  increasing 
function  of  J cs  argument,  and  observing  the  rather  narrow  range  for 
which  2  differs  appreciably  from  zero,  we  conclude  that  ^^(t) 
depends  on  only  a  small  part  of  the  total  creep  compliance  curve 
(e.g.,  Fig.  4).  Since  L  is  the  current  value  of  logarithmic  time, 

Eq.  (65)  shows  that  the  influencial  part  of  the  creep  compliance 
is  in  the  (approximately)  1.2  decade  range  from  (log  t  -  1.2)  to 
log  t.  We  assume  that  over  this  ]  2  decade  range  the  creep  compliance 
can  be  represented  by  the  power  law 

C  (t)  =  cntn  . 
v  1 


(66) 
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Somewhat  arbitrarily,  we  select  the  exponent  n  to  be  the  log-log 


slope  of  the  creep  compliance  at  the  point  L  -  0.48,  which  is  near  the 


centroid  of  the  weighting  function.  Inasmuch  as 


L  -  0.48  ■  log  t  -  log  3.0  ■  log  (t/3.0),  we  see  n  is  the  slope 


of  a  line  drawn  tangent  to  a  double  logarithmic  plot  of  Cv(t)  at 


t/3.  The  coefficient  is  the  value  of  compliance  where  this  tangent 


line  intercepts  the  log  t  «  0  .axis. 


The  actual  value  of  effective  compliance  is  now  readily  derived 


by  substituting  power  law  Eq.  (66)  into  Eq.  (63),  and  treating  n  and 


as  constants  in  performing  the  integration;  there  results 


Ce£Ct)  =1,(1. 


or,  equivalently. 


C  ,(t)  =  C  (A1/n  t) 
ef  v  n 


where 


3/n  T(n  +  1) 


n  4(n  +  |)r(n+f) 


.  n-1  -t  , 
t  e  dt 


is  the  Gamma  function  with  argument  n.  The  factor  A^  is  plotted  in 


Fig.  6  for  the  range  of  n  exhibited  by  polymers  (i.e.,  0  <  n  <  1) ; 


for  most  polymers  0  <  n  <  0.5,  with  filled  polymers  and  plastics 
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usually  having  n  values  at  the  low  end  of  this  range i  Consequently, 
the  effective  compliance  is  commonly  very  close  to  the  creep  compliance 
itself. 

Returning  to  the  displacement  Eq.  (62a),  we  see  that  it  can  be 
written  in  the  simple  form 


v  « - — -  S3/2  A  C  At"  H(C) 

.  r  n  1 

3  TT  rCl 

and  using  Eqs.(62b)  and  (66), 

2  I9  1/0 

v  “ - -  C3/2  C  (t)  H(C) 

3  7r  ya 


where 


t  =  A1/n  C/a 
n 


(71) 


(72a) 


(72b) 


Equation  (72)  is  identical  to  the  elastic  displacement  except 
Cv(t)  appears  in  place  of  Cg. 

This  simple  result  is  easily  generalized  to  the  case  in  which 
higher  order  terms  are  retained  in  the  displacement  expansion. 
Namely,  suppose  the  elastic  solution  is  represented  by  the  sum. 


v  =  l  v  H(C) 
r 


(73) 


in  which 


v  =  A  Cr 
r  r 


C 

e 


(74) 


where  r  >  1  and  A 

r 


is  essentially  constant  during  the  period 


Cl> 
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required  for  the  crack  to  propagate  a  distance  o.  The  viscoelastic 
solution  to  each  term  for  >  t  >  t^  can  be  written  in  the  form 

v  =  A  5r  C™  (At)  (75) 

r  r  et 

where  the  r*"*1  effective  compliance  is 

■t 

C^}(t)  =  r  t"r  Cy(t  -  p)  p*"1  dp  .  (76) 

o 

When  this  compliance  is  expressed  using  a  logarithmic  time  scale 
through  definitions  in  Eq.  (64)  we  obtain  a  result  which  is  identical 
to  Eq.  (65)  except  the  weighting  function  w ,  say,  must  be  changed  using 
the  substitution, 

10£[1  -  10£]1/2  +  |r  10£[1  -  loV'1  5  wr  (77) 

The  weighting  function  for  r  -  |  is  shown  in  Fig.  5,  and  is  seen  to  be 

3 

close  to  that  for  r  =  Indeed,  this  function  is  relatively  insensitive 
to  r;  for  example,  the  maximum  occurs  at 

1  -  -  log  r  (78a) 

with  the  value 

wr(max)  -  j  (1  ”  ~)r”1  (78b) 

am’ 

2 

w^Onax)  -*■  =  0.25  as  r  -*■  00  .  (78c) 

The  general  shape  of  each  wr  is  the  same  as  those  in  Fig.  5.  It  is 
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important  to  realize  that  even  for  r  quite  large  the  essential  part 
of  the  weighting  funccion  lies  in  the  same  £-range  as  the  curves  in 
Fig.  5;  for  example,  with  r  =  10  the  maximum  is  at  t  ■  -1  and  has 
the  value  w^^max)  =  0.26. 

Thus,  the  effective  compliance  can  be  evaluated  quite  accurately 
for  a  wide  range  of  r  using  Eq.  (66)  with  constants  found  from  the 
tangent  to  the  log-log  plot  of  C^(t)  at  t/3.  When  this  power  law 
is  substituted  into  Eq.  (76)  we  find 


C(^(t)  «  X  C  tn 
ef  nr  1 


(79) 


where 


,  .  r  r  (r)  r  (n  +  1) 

nr  (n  +  r)  T(n  +  r) 


(80) 


The  coefficient  X  reduces  to  £q.  (69)  when  we  set  r  =  -r-.  This 
nr  in/  2 

compliance  coefficient  is  drawn  in  Fig.  6  for  four  different  values  of 
r;  X  for  r  =  1  is  shown  because  the  exponent  is  a  lower  bound, 
although  values  of  r  less  than  y  a^e  not  expected  to  appear  in  a 
series  expansion  of  the  exact  displacement  in  view  of  the  analysis 
in  Section  2. 

The  relative  insensitivity  of  Xnj_  to  r,  especially  for  the  usual 
range  of  n  for  polymers  (0  <  n  <  0.5),  leads  immediately  to  a 
simple  solution  to  the  viscoelastic  displacement.  For  example, 
suppose  the  elastic  displacement,  Eq.  (27),  can  be  approximated  by 


a  three-term  series 
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1/2  2  3 

V  “  C  £A3/2C  +  A5/2C  +  A7/2C  1  Ce  H(C)  *  (81a) 

Then  according  to  Eqs.  (75)  and  (79)  an  approximate  viscoelastic 
solution  is 


v  =  51/2[A3/2S  +  A5/252  +  A?/2C3]  Xn  ClAtn  H(C)  (81b) 

and  using  Eqs. (62b),  (66),  and  (72b), 

1/2  2  3  - 

V  =  r  [a 3/25  +  A5/2r  +  A7/2ri  cv(t)  hu>  .  (sic) 


We  have  replaced  the  individual  coefficients  X  with  X  (corresponding 

nr  n 

3 

to  r  *  — )  on  the  assumption  that  the  first  term  in  the  series  is  the 
largest,  and  therefore  X^  i--  more  appropriate  than  a  simple  average 
of  the  three  coefficients.  We  now  replace  the  series  in  Eq.  (8ic)  by 
the  original  expression,  Eq.(27),  and  obtain  the  important  result. 


C  (t) 

V 

VST7-H«) 


Of(C')  2(|r)1/2  -  In 


/ell  '/k 


where  t  is  given  by  Eq.  (72b). 

Although  Eq.  (82)  is  considerably  simpler  than  the  exact 
displacement,  Eq.  (41),  it  is  believed  that  it  is  sufficiently 
accurate  for  our  purposes  regardless  of  the  nature  of  the  failure 
stress  distribution.  We  intend  to  use  this  displacement  to  calculate 
only  the  work  done  on  the  total  failure  zone,  Eq.  (57);  this 
computation  is,  of  course,  an  averaging  process  and  requires  the 
displacement  just  in  the  range  of  0  <  (  <  a.  Moreover,  by  assuming 


39 


several  different  realistic  functions  for  it  is  shown  in 
Appendix  C  that  even  the  one-term  representation,  Eq.  (72a),  is  quite 
good  for  this  purpose. 

In  further  support  of  the  accuracy  of  Eq.  (82),  let  us  consider 

a  few  simple  examples  of  displacement  predictions.  First,  suppose 

oc  is  independent  of  £;  viz.,  o £  =  a  .  The  exact  elastic  displacement 
t  i  m 

is  found  from  Eq.  (26)  to  be 


v  -  j  [iJn  -  (1  -  n)  In  —  ]  HU) 

1  -  ^ 


(83a) 


where  n  =  t/o  and 


A  =  C  o  cx/tt 
e  m 


(83b) 


Expanding  the  logevlthm  in  a  power  series  yields  for  0  <  n  <  1: 

o*  —1/2.1  .121314  .  .  .  ,Q ,  . 

v  =  2A  n  t3'n  +  i5r)  +15n  +  63n  +  "‘l  H(£)  (83c) 

As  a  measure  of  the  error  in  Eq.  (83c)  if  only  three  terms  are 
retained  (as  in  Eq.  (81a)),  we  shall  calculate  the  ratio  of  this 
displacement,  denoted  by  v^,  to  the  exact  displacement,  Eq.  (83a), 
at  n  *  1: 


v3(l) 


=  0.86 


(83d) 


If  only  the  first  term  is  retained  in  Eq.  (83c),  denoted  by  v^,  then 


vx(l) 


(83e) 


40 


As  a  further  comparison  we  calculate  the  ratio  v^(l)/v(l)  given  the 
two  functions 


(i)  of  -  am[l  -  n2] 

(ii)  of  -  a  [1  -  n] 

t  m 


(84a) 


Equation  (38)  may  be  used  to  find  the  one-term  approximation, 
v^(l) ,  with  Eq.  (26)  still  used  to  find  the  exact  displacement  v(l); 
the  result  is 


(i)  Vl(l)/v(l)  -  0.80 

(ii)  Vl(l)/v(l)  -  1.33 


(84b) 


The  first  result,  Eq.  (83d),  although  representing  only  a  14%  error 

in  displacement  at  the  end  of  the  failure  zone,  could  be 

reduced  even  further  by  evaluating  the  three  coefficients  differently; 

e.g.,  a  least  squares  method  could  be  used.  The  above 

results  (and  Appendix  C)  for  a  one-term  representation  give  evidence 

that  the  essential  features  of  the  displacement  at  the  boundary  of 

3/2 

the  failure  zone  are  embodied  in  the  function  4 

Calculation  of  Work  Done  on  the  Failure  Zone:  Equation  (57)  will 
be  evaluated  using  displacement  Eq.  (82).  First,  we  define  an 


v 


(85) 
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where  v  Is  the  displacement  In  Eq.  (82)  and 
—  1  /  n.  • 

t  =  A*'  a/a  .  (86) 

Also,  by  definition,  the  coefficient  C,  and  exponent  n  used  in 
Cy(ta)  and  Eq.  (86)  are  to  be  obtained  from  the  tangent  to  the  double 
logarithmic  plot  of  the  creep  compliance  C^(t)  at  the  time  a/3a; 
note  that  a/a  is  the  time  the  crack  tip  takes  to  propagate  the 
distance  a. 

Now,  divide  the  work  .v.iput  into  two  parts  by  adding  and 
subtracting  v^;  thus 


r  -  r  +  r. 

a  b 


(87a) 


where 


r  = 


3v 

°f «>  ar  d{ 


o 

ra 


af  (0  (v  -  va)  d5 


(87b) 


(87c) 


It  will  be  shown  that  the  second  integral  is  relatively  small  while 
the  first  one  reduces  to  a  simple  function  of  the  stress  intensity 
factor. 

The  Integral  Use  of  Eq.  (85)  yields 


°£<«>  k 


vo 


C  (t) 

I-~-  -  U  v 

C  (t  )  C 

V  O 


dC  • 


(88a) 


I 


0  at  £  =0  and 


Integrate-by- parts  and  use  the  fact  that  * 

C  (t  )  «  C  (t)  at  C  «  o: 
v  a  v 


9a,  C  (t) 

jr1 1-—  - 1)  »„  d5 

5  W 


(88b) 


We  see  that  vanishes  identically  for  two  important  limit  cases: 

(i)  constant  failure  stress  distribution,  a,  =  a  ;  and  (ii)  elastic 

1  m 

continuum.  It  should  be  added  that  in  both  of  these  cases  the  material 

in  the  failure  zone  may  be  highly  nonlinear  and  viscoelastic  without 

affecting  the  vanishing  of  For  cases  other  than  the  above  two 

we  return  to  Eq.  (88a)  and  note  that  the  term  in  curly  brackets  is 

small  in  magnitude  when  v&  is  large  (i.e.,  when  £  ~  at)  and  vice-versa 

(i.e.,  when  £  =  0)  .  Figure  7  shows  the  entire  term  in  curly  brackets 

(normalized  with  respect  to  the  displacement  at  the  end  of  the  failure 
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zone,  v(a))  when  the  displacement  variation  is  given  by  £  ;  this 

use  of  a  one-term  representation  of  displacement  is  not  crucial  to 
the  argument  as  other  cusp-type  displacement  variations 
would  lead  to  the  same  general  behavior  shown  in  Fig.  7. 

We  have  used  the  same  constants  (n,C^)  in  both  compliances  in 
Eq.  (88a),  which  yields 


rCv(t)  ,,  n,  3/2  ,  N 

[ : 1]  v  =  -  [l  -  n  ]  n  v(a) 

c  (t  )  a 

v  a 


(89a) 


Now,  r,  is  the  work  done  by  a.  acting  through  a  displacement  equal  to 
b  I 
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the  entire  quantity  in  Eq.  (89a).  By  comparison,  Eq.  (87b)  shows 
that  Ta  is  the  work  done  by  acting  through  the  displacement  v^, 

3/2 

va  =  n  v(a)  .  (89b) 

Examination  of  the  displacements  in  Fig.  7  reveals  clearly  that 

I*  >>  Tb  unless  acts  mainly  in  the  neighborhood  of  n  ■  0.  Such 

behavior  of  cannot  exist  in  view  of  the  definition  of  the  failure 

zone;  viz.,  a  is  defined  as  the  distance  over  which  only  significant 

values  of  o^.  act,  and  if  acted  only  at  the  far  left  of  Fig.  7, 

a  would  be  reduced  accordingly.  However,  it  is  of  interest  to  point 

1/2 

out  that  for  even  the  rather  rapidly  decaying  stress  o^  =  0^(1  “  h  ) 

the  ratio  T  /T  is  n/ (n  +  2)  and  this  ratio  is  less  when  the  actual 
b  a 

3/2 

displacement  variation  is  used  in  the  calculation  instead  of  n 

(If  one  wished  to  retain  T^,  its  evaluation  using  even  a  crude  estimate 

of  probably  would  be  adequate  considering  its  smallness.) 

As  one  final  point  concerning  T^,  we  should  mention  that 
quantities  n  and  generally  will  vary  with  n,  rather  than  being 

constant  as  assumed  in  deriving  Eq.  (89a).  However,  by  changing 
Eq.  (88b)  to  a  logarithmic  scale,  log  n,  and  using  arguments  analogous 

to  those  used  to  derive  Eq.  (67),  it  can  be  shown  that  is  essentially 

-?  _  .... 

independent  of  t  except  for  the  range  10  t  <  t  <  t^;  this 
observation  together  with  the  typically  weak  dependence  of  and  n 
on  log  t  justifies  our  use  of  Eq.  (89a). 

The  Integral  T  :  We  shall  write  out  r  ,  Eq.  (87b),  by  using 
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Eq.  (82)  and  definition  Eq.  (85);  after  performing  the  differentiation 
there  results,  finally. 


r 

a 


C  (t  ) 
v  a 

-a 

'a  0f(O  Of(C) 

2  7T 

0 

n  G  JV 

0 

Z  -  v 

»  j 

di'  d£ 


(90a) 


Note  that 


5  +  C 


5  -  5 


-) 


(90b) 


Substitution  of  Eq.  (90b)  into  (90a)  and  use  of  Eq.  (18)  for  the 
stress  intensity  factor  yields 


TT  ■*  2 

7  c  (t  )  r  +  z 

4  v  a  o 


(90c) 


where 


Z 


c  (t  ) 

v  a 

■a 

of(5)  of(£') 

'  ' 

t 

4  tr 

0 

.  * 

6  -  5’ 

*  > 

(90d) 


The  integral  Z  is  identically  zero.  This  result  is  clearly  seen  by 
referring  to  Fig.  8  and  observing  that  the  integrand,  J?(£,  ,£ ') ,  say, 
is  skew- symmetric  with  respect  to  the  line  §  =  £; ' ;  viz.,  by  direct 
substitution  one  finds 


=  -  4tt2,ep 


(90e) 


Governing  Equations  for  Crack  Growth:  Neglecting  on  the  basis 
of  the  above  arguments  and  bringing  together  Eqs.  (87a)  and  (90c) 


i 


yields  the  principal  result  of  this  paper: 


C  (t  ) 
v'  a 


This  equation  is  identical  to  that  for  the  critical  stress  intensity 

factor  of  an  elastic  material  (see  Eq.  (59))  except  the  linear 

viscoelastic  creep  compliance  Cv(t^)  appears  in  place  of  the  elastic 
2 

constant,  4(1  -  v  )/E.  For  completeness  we  record  here  the  definition 


of  t  : 
a 


t  =  X1/n  a/a 
a  n 


where  a  is  found  frcm  Eq.  (22), 


2  »2 

it  N 

_ o 

2  2 

a  I 

m  1 


and  A  '  is  obtained  from  Fig.  6,  with  the  value  of  n 


(92a) 


(92b) 


simply  related  to  creep  compliance,* 


d  log  Cv(t) 


d  log  t 


(92c) 


at  the  time  t  =  a/3a.  Experimental  and  theoretical  determination 

of  fracture  properties  T  and  a  I  ,  and  their  possible  velocity 

m  1 

dependence,  will  be  discussed  a  little  later;  for  the  time  being  we 
shall  treat  them  as  known  quantities. 

The  general  character  of  continuous  crack  propagation  in 
viscoelastic  media  can  now  be  described.  First,  we  note  that 
Cy(ta)  is  a  continuous,  monotonic,  increasing  function  of  t  having 

*In  most  problems,  the  value  A  ^  =1/3  obviously  will  be  acceptable. 


the  range  of  variation, 


Cv(0)  <  CvCta)  <  Cv(-)  .  (93) 

Thus,  Eq.  (91)  admits  a  single-valued  continuous  solution  tQ,  given 
Nq  and  f,  for  the  following  rang:;  of  stress  intensity  factors 


N  <  N  <  N 
oe  o  og 


(94a) 


where 


^  /  4  r 

oe  “  y  TT  Cv(°°) 


(94b) 


and 


"og 


4  r 


c  (0) 

v 


(94c; 


These  limiting  stress-intensity  factors  are  identical  to  so-called 


critical  stress-intensity  factors  of  elastic  media  having  elastic 


2  2 

constants  4(1  -  v  (#o))/E(00)  and  4(1  -  v  (0))/E(0),  respectively; 


the  relation  between  C  (°°) ,  C  (0)  and  these  elastic  constants 

v  v 


follows  from  application  of  the  initial  and  final  value  theorems  of 

Laplace  transform  theory  [17]  to  Eq.  (45). 

Henceforth,  N  and  N  will  be  called  the  equilibrium  and 
oe  og  — 1 - 

glassy  factors,  respectively;  the  equilibrium  factor  for  a  non- 


crosslinked  polymer  is  zero  since  -  °°. 


Returning  to  Eqs.  (91)  and  (92),  we  see  the  behavior  described 
above  implies  the  instantaneous  crack  tip  velocity  a  is  a 
monotonically  increasing  function  of  N  in  the  range  in  Eq.  (94); 
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although  T  and  1^  may  depend  on  velocity  it  is  doubtful  this 

dependence  would  alter  this  general  behavior.  As  Nq  approaches  N 

from  below,  a  becomes  unbounded,  while  as  N  approaches  N  from  above 

o  rr  oe 

• 

a  approaches  zero.  It  is  concluded,  therefore,  that  a  crack  will  not 

propagate  if  Nq  <  N  and  will  propagate  at  high  velocities  (limited 

of  course  by  effects  of  wave  action  which  have  not  been  taken  into 

account  in  this  paper)  when  N  >  N 

o  -  og 

Another  point  of  interest  is  that  instantaneous  tip  velocity 

depends  on  the  current  value  of  N  but  not  on  past  values.  Recall 

o 

also  that  a  is  actually  the  velocity  normal  to  the  curve  defining 
the  crack  tip,  and  consequently  Eq.  (91)  enables  one  to  compute  the 
time-dependent  shape  and  size  of  a  crack  in  terms  of  local  stress 
intensity  factors. 

In  order  to  solve  Eq.  (91)  explicitly  for  tip  velocity  it  may  be 
desirable  to  reintroduce  the  power-law  form  C^(t)  =  C^tn;  Eqs.  (91) 
and  (92)  then  yield,  with  a  =  da/dt: 


da 

dt 


C  A  7r(2n+l)  1/n 
1  n 
4  r 


N2(l+l/n) 

o 


2  2 
o 

m  1 


(95) 


(In  the  development  of  the  theory,  a  was  defined  to  be  the  velocity 
at  the  time,  t^,  which  is  the  time  at  which  the  crack  reaches  a 
generic  location  x;  for  notational  simplicity  we  now  drop  the 
subscript  "1"  on  the  time  variable.)  In  general,  this  is  a  first  order 
nonlinear  differential  equation  for  crack  length.  If  n  is  not 
constant,  it  has  to  be  expressed  in  terms  of  N  and  a  through 
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Eqs .  (92b)  and  (92c).  (In  solving  Eq.  (95)  numerically,  a  considerable 
reduction  in  computational  time  may  be  achieved  by  letting  n  be 
piecewise  constant  over  approximately  one-decade  intervals  in  a/a, 
depending  on  the  amount  of  curvature  in  the  log-log  plot  of  creep 
compliance;  this  is  especially  true  when  Eq.  (95)  adtu.us  an  analytical, 
closed-form  solution  with  n  constant.) 

The  solution  of  Eq.  (91)  may  be  facilitated  by  incorporating 
the  initial  value  of  compliance,  Cv(0),  ir.  the  analytical  representa¬ 
tion  of  Cv(t);  we  do  this  by  means  of  the  generalized  power  law; 

Cv(t)  =  Cq  +  C2tm  (96) 

where  =  Cv(0);  m  and  C 2  are,  respectively,  the  slope  and  the 
log  t  =  0  intercept  of  a  line  drawn  tangent  to  a  double-loga'  xthmic 

plot  of  [C  (t)  -  C  ]  at  time  t.  The  main  advantage  of  this 

vo 

representation  over  simple  power  law  Eq.  (66)  is  that  m  and  C 2 
are  constant  over  many  decades  of  time  for  numerous  plastics  and 
elastomers  [  1  ];  for  example,  the  above  form  of  power  law  fits  the 
uniaxial  data  in  Fig.  4  extremely  well  for  -  00  <  log  t  <  -2;  in  fact, 
except  for  a  small  portion  of  the  curve  close  to  log  t  =  -2,  the 
entire  compliance  can  be  represented  by  the  generalized  power  law  up 
to  its  intersection  with  D(°°)  and  by  the  constant  D(°°)  at  greater 
times.  (Observe  that  the  condition  of  small  curvature  used  in 
developing  the  theory  (e.g.,  Eq.  (67)),  which  requires  the  slope  n 
to  be  essentially  constant  over  approximately  1.2-decade  intervals,  is 
not  particularly  well  satisfied  close  to  log  t  =  2.  It  is  important 
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to  realize  that  this  curvature  will  produce  some  error  in  Eq.  (91) 
only  when  log  (a/3a)  =  2,  but  will  not  cause  error  at  other 
velocities . ) 

One  can  easily  show  that  if  generalized  power  law  Eq.  (96)  had 
been  used  in  the  development  of  the  theory,  instead  of  Eq.  (66), 

Eq.  (91)  still  would  have  been  obtained.  However,  one  must  replace 
n  by  m  in  Eqs.  (92a)  and  (92c)  and  in  Fig.  6;  also,  replace 
Cv(t)  by  [Cv(t)  -  C^]  in  Eq.  (92c).  We  now  solve  Eq.  (91)  for  tip 
velocity  after  drawing  upon  Eq.  (96)  and  definition  Eq.  (94c): 


da 

dt 


Cl  Z2"*1) 

z  m _ 

No 

[4F[1  --2-] 

\i  ™ 


1/m 


N2 (1+1/m) 

o 


(97) 


2  2  2  2 

for  N  <  N  .  We  see  that  when  N  /N  <<  1  this  result  reduces 
o  og  o  og 

to  Eq.  (95);  for  this  range  we  may  set  C ^  and  ra  =  n.  Equation  (97) 

also  shows  clearly  that  the  velocity  becomes  unbounded  when 

approaches  N  from  below, 
og 

Fracture  Properties:  The  basic  Eq.  (91)  was  derived  without 
explicitly  imposing  any  significant  restrictions  on  the  nature  of 
the  failure  stress,  o^.  Now  we  shall  deduce  some  special,  but  still 
relatively  general ,  representations  in  order  to  be  able  to  make 
explicit  predictions  of  crack  growth  and/or  determine  pertinent 
fracture  properties  from  experimental  data. 

First  it  should  be  pointed  out  that  the  assumed  geometry  of 
the  crack  cross-section.  Fig.  1,  probably  limits  us  to  functions 
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o^(£)  which  are  either  independent  of  £  or  are  decreasing  functions 

of  £.  For  if  a increased  appreciably  with  £  (such  as  for  a  metal 

with  significant  strain  hardening),  the  continuum  adjacent  to  the 

failure  zone  would  be  subjected  to  stresses  outside  of  its  range  of 

linear  viscoelasticity,  thereby  invalidating  a  basic  assumption  of 

the  theory.  This  limitation  is  not  believed  to  be  a  serious  one 

for  most  filled  and  unfilled  polymers;  certainly  the  theory  will 

be  suspect  for  some  materials,  such  as  natural  rubber  in  which 

stress-  and  rate-dependent  crystallization  forms  a  strong,  anisotropic 

crack  tip  reinforcement  [  2  ] .  On  this  basis  it  will  be  assumed 

that  the  maximum  value  of  a  ,  which  is  denoted  by  a  ,  occurs  at 

i  m 

the  tip  or  at  least  very  close  to  it. 

Let  us  now  determine  the  principal  geometric  and  mechanical 
variables  affecting  the  dimensionless  failure  stress  distribution 
f (-  at  a  generic  location,  £.  In  general,  this  stress  will 

depend  on  the  displacement  v  at  (  as  well  as  the  t*._.  variation 
of  displacement  as  it  increases  from  zero  to  the  maximum  value  v 

m 

at  £  =  a.  This  time  variation  of  displacement  is  defined  by  the 
tip  velocity  a  (which  already  has  been  assumed  constant  during 
the  time  an  element  of  material  is  stretched  from  its  initial  length 
at  £  =  0  to  complete  failure  at  a),  the  shape  of  the  interface 
jetween  the  linear  continuum  and  failure  zone,  and  the  maximum 
displacement  v  ;  we  shall  express  the  interface  shape  using  the 
single  parameter  s g/ ,  which  is  intended  to  represent  the  set  of  all 
values  of  v/v  over  0  <  £  <  a.  The  rate  at  which  the  tensile  stress 
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in  the  x-direction  within  the  failure  zone  decreases  from  its  value 

of  o  (o)  (=  o  (o))  at  C  =  0  to  zero  could  affect  the  rate  of 
x  y 

degradation  and,  therefore,  o^;  e.g.,  the  triaxial  tensile  stress 

state  near  the  tip  may  facilitate  the  production  of  cavities  in  the 

failure  zone  for  some  materials.  It  will  be  assumed  that  this 

effect  is  determined  by  the  overall  failure  cone  dimensions,  v  and 

m 

a,  tip  velocity,  a,  and  interface  shape  of. 

One  might  also  argue  that  the  crack  tip  thickness  d^  (see  Fig.  3) 
is  not  a  constant,  and  therefore  should  be  explicitly  included  in 
the  set  of  parameters  affecting  f.  However,  implicit  in  the 
assumption  of  linear  viscoelastic  behavior  for  the  material  outside 
of  the  failure  zone  is  that  the  strains  are  small,  and  this  restriction 
includes  the  strain  at  the  tip,  5=0;  within  this  context  d^  is 
essentially  constant  and  no  further  mention  of  this  thickness  need 
be  made. 

On  the  basis  of  the  above  remarks  and  dimensional  considerations, 
the  dimensionless  distribution  f  can  be  written  as  follows: 


““  “/V"  the  Influence  ef  th.  trlaxlil 

streee  etate  on  the  degradation  proceea. 

°tdCr  C°  eSt8bll8h  variables  affecting  th. 
streas  consider  a  nateriel  point  along  the  line  of  crack 
prolongation.  As  the  crack  tip  approaches,  the  stresses  (o  ,  o  „ 
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-tly  speaking  o,  depends  on  stress  history  as  shown  i„  (98b) 

"  CO"Blder1"8  -  ^  insensitivity  of  Poisson.,  ratio 
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Xt  18  of  interest  to  mention  that  th* 

“  the  above  representation  of 

f  °a  F  contains  as  a  special  case  the  fracture  theory  of 

Barenblatt  et  al.ti8J>  „hlch  ls  ba8ed  „„  ^  ^  ^ 

ory  to  the  failure  zone  in  an  otherwise  linearly  elastic 
material.  However,  we  shall  «  their  special  form  of 


as  it 
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is  derived  from  a  rate  equation  for  bond  failure  which  is  probably 
over-simplified  [2]  and  also  predicts  the  physically  incorrect 
result  of  a  non-zero  tip  velocity  for  a  vanishingly  small  stress 
intensity  factor;  while,  such  behavior  could  exist  in  noncrosslinked 
polymers  above  their  glass-transition  temperature,  they  do  not  satisfy 
the  elasticity  assumption  made  by  Barenblatt  et  al.[18]  when  a  0. 

Returning  to  the  discussion  leading  up  to  Eqs.  (98a)  and  (98c), 
we  may  conclude  that  the  fracture  energy  (which  is  the  work  done  by 


a  acting  through  the  displacement  v  )  depends,  at  most,  on  the 
i  m  ' 

nondimensiunal  parameters  shown  in  Eq.  (98a)  except  for  v/v  ;  viz., 

m 


m 


(99) 


The  quantities  T,  o  ,  I  ,  and  F,  by  themselves  or  in  combination, 

m  1 

will  be  called  fracture  properties.  For  the  sake  of  argument,  let 
us  now  suppose  that  these  fracture  properties,  as  well  as  Cv(t),  are 
given  functions  of  their  arguments.  Equations  (91),  (92),  and 
Eq.  (82)  provide  a  set  of  five  equations  from  which  the  five 
quantities  n,  t  ,  a,  a,  and  v  can  be  found,  in  principle,  given 
N^;  note  that  the  shape,  js/*,  is  defined  by  the  variation  of  v  with 
i  through  Eq.  (82).  It  is  therefore  concluded  that  the  instantaneous 
crack  tip  velocity  depends  on  only  the  instantaneous  stress  intensity 
factor;  neither  the  history  of  this  factor  nor  the  history  of  the 
stresses  in  the  continuum  enter. 


This  conclusion  has  an  important  practical  implication.  Suppose 
that  the  function  Nq  =  h(a)  has  been  obtained  experimentally  over  the 
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velocity  range  of  interest.  Then  to  predict  crack  growth  due  to  a 
stress  intensity  factor  which  has  different  time-dependence  than  used 
in  the  experimental  determination  of  h(a),  one  does  not  have  to 
introduce  any  new  functions  or  fracture  theory  per  se.  Rather,  all 

one  need  do  is  integrate  the  first  order  differential  equation 

•  • 

Nq  *  h(a),  where  may  possibly  depend  on  a  as  well  as  other 

time-dependent  parameters.  Of  course,  use  of  the  previously  derived 

theoretical  relationships  may  greatly  facilitate  the  evaluation  of 

h(a)  from  experimental  data  (especially  if  T  and  a  I  are  constant) 

m  1 

and  simplify  the  determination  of  its  dependence  on  such  important 
factors  as  chemical  aging  and  temperature;  incorporation  of  these 
factors  will  be  discussed  shortly.  For  the  three  materials  studied 
later  in  Section  4  it  is  found  that  T  and  o  I,  can  indeed  be 
taken  as  constants;  this  result  will  be  seen  to  provide  a  simple 
analytical  representation  of  h(a)  in  terms  of  creep  compliance 
and  enable  some  rather  useful  conclusions  to  be  drawn  concerning 
crack  growth  and  failure  under  realistic  loading  conditions. 

Finally,  it  should  be  clear  that  if  only  the  function  h(a)  is 
known,  it  is  not  possible  to  establish  the  separate  fracture 
properties  appearing  in  the  theory.  For  example,  we  could  arbitrarily 
assume  F  =  1  (or  1^  “  2)  and  select  any  constant  value  for  and 
still  adjust  the  theory,  Eq.  (91),  so  as  to  fit  the  function  f(a); 
this  fitting  would  be  accomplished  by  absorbing  as  much  velocity 
dependence  in  T  as  necessary  by  setting  Nq  ■  h  and  combining  Eqs.  (91) 
and  (92)  so  as  to  find  T  from  the  equation. 
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We  are  thus  led  to  an  important  conclusion  concerning  fracture 

experiments.  If  it  is  desired  to  use  these  experiments  to  establish 

the  rate  dependence  of  T  and  a  I,,  if  indeed  any  exists,  one  must 

m  1 

measure  both  crack  velocity  and  failure  zone  length  as  functions  of 

N  . 
o 

Effect  of  Environmental  Changes  and  Aging;  The  creep  compliance 
and  fracture  properties  T  and  which  appear  in  the  basic 

fracture  Eqs.  (91)  and  (92)  can  bs  expected,  in  general,  to  vary 
with  the  environment  through,  for  example,  dependence  on  temperature 
and  relative  humidity.  These  properties  may  also  vary  directly  with 
age  due  to  such  phenomena  as  post  cure  reactions  and  oxidation.  As 
long  as  the  resulting  changes  in  Cv»  T,  and  are  relatively  small 

in  the  generic  time  interval  a/a,  we  can  incorporate  this  time 
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variation  in  properties  directly  into  Eqs.  (91)  and  (92),  as  well 
as  Eqs.  (95)  and  (97),  without  having  to  alter  the  form  of  the 
equations;  integration  of  the  appropriate  equation,  taking  into 
account  the  time-varying  properties,  will  yield  instantaneous  crack 
size  as  a  function  of  environment  and  aging  histories. 

As  a  point  of  caution  it  should  be  recalled  that  the  crack 
propagation  theory  was  developed  under  the  assumption  that  Nq  is 
independent  of  viscoelastic  properties  of  the  continuum.  With 
sufficiently  rapid  changes  in  the  environment  and/or  spacially 
nonuniform  aging,  the  resulting  inhomogeneity  may  give  rise  to 
significant  dependence  of  upon  the  viscoelastic  properties. 
Consequently,  the  theory  will  not  be  valid  during  those  periods  for 
which  such  inhomogeneities  exist. 

A  study  of  the  possible  ways  in  which  the  fracture  properties 
may  change  is  beyond  the  scope  of  this  paper.  But  we  will  comment 
briefly  on  temperature  dependence  of  the  creep  compliance  as  a  large 
amount  of  experimental  data  exists  on  this  property.  For  many 
crystalline  and  amorphous  polymers,  below  and  above  the 
glass-transition  temperature,  T  ,  temperature  dependence  of  the 
compliance  is  characterized  under  isothermal  conditions  by  the 
simple  equation  [  1  ] , 

Cy(t,  T)  =  Cq  +  AC(£)/aG  (101) 


where 


Cq  =  Cq(T)  ■  initial  compliance 
4C«)  =  (Cv  -  co)  *e 
5  =  t/aT  -  reduced  time 
aG  =  "  compliance  factor 

aT  "  ax^T^  *  time-scale  factor 

Thus,  three  functions  of  temperature  (C  ,  a  ,  a„)  plus  one  function 

of  reduced-time  (AD(£))  define  the  compliance  C(t,  T) .  It  should 

be  added  that  most  of  the  data  used  to  establish  the  form  of 

Eq.  (101)  are  from  uniaj.ial  tensile  tests;  however,  this  equation 

should  be  quite  accurate  for  C  as  well  since  the  uniaxial  creep 

2 

compliance  and  differ  by  only  the  factor  (1  -  v  ) • 

Let  us  assume  the  common  power  law  form  fcr  AC(£)» 

AC(£)  -  CRCm  (102a! 

where  CR  and  m  are  constants.  Upon  comparing  Eqs.  (96)  and 
(102a)  we  see  that  is  a  function  of  temperature;  viz., 

C2  -  Cg/a™  .  (102b 

This  coefficient,  as  well  as  Cq(T),  can  be  substituted  directly 
into  Eq.  (97)  in  order  to  predict  crack  growth  under  different 
constant  or  transient  temperatures. 

When  amorphous  polymers  are  in  their  glass-to-rubber  transition 
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:  •  one  can  neglect  the  temperature  dependence  of  and  Cq, 

ur.  aerefore  write 

Cv  5  Cv(5)  (103) 

where  £  =  t/a^.  Any  material  whose  temperature  dependence  enters 
entirely  through  the  time  scale,  as  in  Eq .  (103),  is  called 
"thermorheologically  simple" [13] . 

Prediction  of  Initial  and  Intermittent  Crack  Growth 

The  preceding  analysis  enables  one  to  predict  crack  growth 
under  a  constant  or  varying  stress  intensity  factor  during  those 
periods  for  which  Nq,  T,  ar.d  do  not  vary  appreciably  over  each 
time  increment  a/a.  Of  course,  even  if  there  are  intervals  during 
which  this  condition  is  not  met,  the , theory  is  still  valid  for 
predicting  growth  in  each  intervening  period  when  is  sufficiently 
large. 

Prediction  of  the  time  at  which  a  crack  starts  to  grow  (i.e., 

the  fracture  initiation  time,  t.)  and  growth  when  N  varies 

i  o 

significantly  during  the  time  increment  a/a  (as  in  low  stress  level 
fatigue)  is  accomplished  by  bringing  together  the  convolution 
representation  of  displacement,  Eq.  (41),  the  finite  stress 
condition,  Eq .  (22),  and  the  fracture  energy  Eq.  (54)  or  Eq.  (55). 
The  analysis  may  be  considerably  more  involved  than  the  preceding 
one  of  continuous  growth  since  simplified  displacement  Eq.  (82) 
cannot  be  used  and  the  energy  T  and  failure  stress  distribution 


... 
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Failure  of  the  material  at  *  -  0  occurs  when  v  reaches  the 
value  v  corresponding  to  the  displacement  for  which  the  energy 
input  is  equal  to  T;  specifically,  from  Eq.  (54), 


F  =  a  v 
m  m 


(107) 


ng  t  t±  and  v  -  vm  in  Eq.  (105)  and  substituting  this  result 
into  Eq .  (107)  we  obtain 


r  =  lNo<ti)  Cv2,<V 


(108) 


which  is  identical  to  the  condition  for  the  onset  of  crack  growth 

elastic  body  except  the  secant  compliance,  Eq.  (106),  replaces 
4(1  -  v  )/E  (see  Eq.  (59)).  Referring  to  Eq.  (94),  we  conclude  that 

Ci  0  if  No(0+)  >  Nog,  and  that  fracture  does  not  initiate  if 

N  <  N  . 
o  oe 

Secant  compliance  Eq.  (106)  can  be  rewritten  as  follows: 


cv2,<t>  '  Cy(t)  -  ~L_ 

N  (t) 

O  J  o 


d  N2(t) 


,Cv(t)  '  cv(t  "  t)]  — -  dt  (1Q< 


oince  Nq  and  are  nondecreasing  functions  of  time,  this  equation 


implies 


c!2,(t)  <  c  ft)  . 
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Hence,  if  the  creep  compliance  itself  were  used  in  the  elastic-like 
Eq.  (108)  instead  of  the  secant  compliance,  the  predicted  initiation 
time,  tci>  say,  would  satisfy  the  condition 


fcci  -  Ci  *  (109c) 

The  equality  results  if  Nq  is  constant  for  t  >  0. 

When  Nq  is  a  power  law  in  time,  the  secant  compliance  can  be 
obtained  from  the  previous  result  for  effective  compliance,  Eq.  (79). 
Namely,  if  Nq  s,  t  ,  where  p  is  a  positive  constant,  we  find 


A  G 

nr  v 


(t) 


(110) 


where  ^  is  given  by  Eq.  (80)  and  r  -  2u.  Figure  6  can  be  used 
to  estimate  A^  if  1/2  <  y  <  7/4. 

The  fracture  initiation  time  predicted  by  Eq.  (108)  was  derived 
assuming  of  *  am  -  constant.  However,  the  results  in  Appendix  C 
provide  evidence  that  Eg.  (108)  is  at  least  approximately  valid  even 
varies  with  ( .  Specifically,  when  approximate  displacement 
Eq.  (49)  is  used  along  with  Eq.  (22)  to  predict  vm,  and  it  is 
assumed  that  the  shape  of  the  failure  rone  is  the  same  as  in  Eq.  (72) 
(via.,  v  v  c"+3/2)  and  that  of  -  »_?(»/»_),  there  results 


r  '‘iiVV  c 


(2) 


(ti> 


(111a) 


where  A^  is  given  by 


A 


i 


f  1 

F(p)  dp 
o 


(111b) 


din 


. . . 


-  -  -  •  -  a., - icwjwi 


The  coefficient  A.  turns  out  to  be  equal  to  the  ratio  T  /T  in 
i  ap 

Eqs.  (04)  and  (C— 5)  when  is  given  by  Eqs.  (C-l)  and  02), 

respectively.  In  view  of  the  results  shown  in  Fig.  C-l,  it  is 

tentatively  concluded  that  Eq.  (108)  can  be  used  to  predict  fracture 

initiation  time  even  when  varies  with  C.  (There  are  cases, 

however,  when  t^  is  very  sensitive  to  T  and  therefore  very  sensitive 

to  the  variation  of  o^.  For  example,  this  situation  exists  when 

(2) 

is  a  constant  within  the  range  in  Eq.  (94a)  and  O  is  a  weak 
function  of  time.) 

Let  us  now  consider  the  problem  of  predicting  fracture  initiation 

time  when  we  remove  the  assumption  that  is  nondecreasing  but,  for 

simplicity,  retain  the  assumption  that  af  =  am>  Equation  (104) 

implies  the  position  of  the  crack  tip  a(=  a)  will  move  to  the 

left  (a  <  0)  upon  unloading  (dN  /dt  <  0)  if  a  does  not  decrease  at 

o  m 

least  in  proportion  to  N^;  result  a  <  0  is  physically  unacceptable 

since  it  would  imply  complete  healing  of  the  material  in  the  failure 

zone  and,  moreover,  would  invalidate  Eq.  (41)  since  it  is  restricted 

to  a  >  0.  At  any  time  (0  <  t  <  t^)  during  the  loading  or  unloading 

© 

periods  the  displacement  is  given  by  Eq.  (41)  as  long  as  a  >  0  and 

the  crack  faces  do  not  contact  one  another  (which  could  occur  if 

N  <  0) .  At  x  =  0  and  with  or  =  a  (t)  and  C  =  a  =  a(t),  we  find 
o  r  m 

ft 

v  =  —  C  (t  -  x)  4“  [°  (T)  <*(t)]  dt  .  (112) 

it  v  di  m 

■  o 

One  physically  reasonable  selection  for  a  (x)  and  a(x)  would  be 

a  *  constant  when  dN  /dt  >  0  and  N  >  N  (max),  where  N  (max)  is  the 
m  o  -  o-o  o 

iargest  N  up  to  the  present  time,  and  a  =  constant  whenever  N  <  N  (max), 
o  oo 
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The  analysis  will  not  be  pursued  further  in  this  paper.  However,  it 
is  to  be  noted  that  the  resultant  work  input  to  the  failure  2one 
at  t  *  t  should  be  calculated  using  Eq.  (55),  rather  than  Eq.  (107), 
since  o  is  time  dependent. 

Initial  Velocity:  Calculation  of  tip  velocity  during  the  first 

stages  of  crack  growth  will  be  illustrated  using  the  assumptions 

that  N  is  nondecreasing  and  a,  =  a  =  constant.  Motion  of  the  crack 
o  f  m 

tip  prior  to  t^  is  obtained  directly  from  Eq.  (104)  since 
da/dt  =  da/dt: 


da 

dt 


dN 
_ o 

dt 


0  <  t  <  t. 


(113) 


m 

Immediately  after  fracture  initiates, the  tip  velocity  becomes 
dependent  on  the  creep  compliance  of  the  continuum.  This  velocity 
can  be  calculated  from  Eq.  (91)  if  Nq  does  not  change  appreciably 
during  the  time  it  takes  for  the  crack  to  propagate  the  distance 
a  in  Eq.  (104).  Inasmuch  as  Nq  satisfies  Eq.  (108)  at  t^,  comparison 
of  Eqs.  (91)  and  (108)  yields  an  implicit  equation  for  a, 

C  (A1/n  a/a)  -  C(2)(t.)  (114) 

v  n  v  l 

where  t.  <  t  <  t.  +  a/a.  If,  for  example,  N  is  constant  for  t  >  0 
i  i  o 

and  N  <  N  <  N  , 
oe  o  og 


(ti>  ■  Cv(tl> 


(115a) 


and  therefore, 
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a 

a 


t  /A1/n 
ti/An 


(115b) 


or 


a 


2  2 

*V(t.) 

O  1 


(115c) 


m 

which  is  the  velocity  immediately  after  the  fracture  initiates  at 

time  t  .  Referring  to  Fig.  6,  it  is  seen  that  0.29  <  A^n  <  0.40 
l  -  n  - 

when  0  <  n  <  1;  thus,  according  to  Eq.  (115b)  the  time  for  the 
crack  to  propagate  the  distance  a  is  approximately  three-times  the 
initiation  time.  In  deriving  Eq.  (115b)  it  was  assumed  that  the 
value  of  n  obtained  from  the  creep  compliance  at  t  =  a/3a 
(see  Eq.  (92c))  is  the  same  as  that  at  t^;  this  is  justified  since 


Eq.  (115b)  and  the  average  value  of  A^n  ^  1/3  imply  a/3a  *  t.  . 

n  l 

Equation  (114)  is  based  on  the  assumption  that  the  fracture 
energy  required  to  produce  the  initial  failure  at  x  =  0  is  equal  to 
that  under  continuous  growth.  Whether  or  not  this  assumption  is 
valid  probably  depends  on  the  material  and  loads  since  the  former 
energy  could  depend  on  loading  history  when  t  <  t^  and  the  latter 
energy  may  be  velocity  dependent.  Furthermore,  an  element  of 
material  in  the  failure  zone  at  x  =  0  is  unstressed  in  the 
x-direction  for  t  >  0,  while  in  the  continuous  growth  case  an  element 
of  material  in  the  failure  zone  is  under  triaxial  stress  at  least 
part  of  the  time. 

Low  Stress  Level  Fatigue:  When  Eq.  (91)  cannot  be  used  to 
predict  tip  velocity,  we  must  return  to  the  convolution 


>  il  ri~  -Tirrai**1^ 
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representation  of  displacement,  Eq.  (41),  in  order  to  establish  a 
criterion  for  crack  growth.  Such  may  be  the  case  in  the  important 

problem  of  cyclic  loading  when  the  maximum  stress  intensity  factor 

\ 

is  large  enough  to  cause  crack  growth  while  the  minimum  factor  is 
so  low  that  either  the  crack  is  stationary  during  a  portion  of  each 
cycle  or  a/a  is  not  small  compared  to  the  vibration  period. 

The  method  of  analysis  used  by  Knauss  and  Dietmann  [7] 
apparently  could  be  applied  to  the  present  theory;  they  expanded 
the  stress  intensity  factor  in  a  power  series  in  time  over  the 
interval  in  which  the  crack  moved  a  distance  a,  where  a  was  constant 
in  their  theory.  Here  it  would  be  necessary  to  account  as  well  for 
the  fact  that  a,  and  possibly  a^,  vary  with  time  over  each  cycle, 

but  such  an  analysis  will  not  be  pursued  in  this  paper. 

4.  Applications 

A  Sheet  Under  Constant  Stress 

Consider  a  centrally  cracked  plate,  such  as  shown  in  Fig.  9. 

Assuming  the  crack  length,  2a,  is  large  compared  to  sheet  thickness 

but  small  compared  to  the  in-plane  dimensions,  the  stress  intensity 
factor  is  [5 ] , 


where  o  is  the  applied  stress.  Crack  growth  and  failure  time  will 
be  predicted  using  Eq.  (97)  under  the  assumption  that  the  creep 
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compliance  parameters,  and  m,  and  ie  fracture  properties,  r  and 
o^I^,  are  constant.  Substitution  of  Eq.  (116)  into  Eq.  (97), 
assuming  o  is  applied  at  t  *  0  and  is  constant  thereafter,  and 
integrating  yields 


,  o  2  I2  r  C  a/n  a/a°  a  . 

t  =  L.  -5—1  .  -2  r 8—  .  iil/m  ±1 

2  2  L  Cj  la  Y  J  Y 

ir  o  v  m  2J  ,  o'  T 

1 


(117a) 


where  a  is  defined  as 

g 


_  sr 

3g  =  - 7 

8  tt  C  a 
o 


(117b) 


When  the  central  crack  reaches  length  2a^  complete  failure  of  the 

sheet  occurs;  this  point  follows  from  Eqs.  (91),  (94c),  and  (116), 

which  imply  a  -*•  «  as  a  -*  a  .  The  time  at  which  a  =  a  is  called 

8  g 

the  failure  time,  t^,  which,  after  changing  the  intergration  variable 
in  Eq.  (117a),  becomes 


o2  I?  t  C  a/m  f1 


9  a  i.  l 
I  m  1  o 

2  2  A  C„ 

tt  a  '■m2-' 


(1  -  u)1'”  u-(1+1/m)du  .  (118) 


a  /  a 
o  g 


It  is  of  interest  to  notice  that  when  a  /a  >>  1  and  0  <  m  <  1, 

go  - 

the  integrand  in  Eq.  (117a)  goes  quickly  to  zero  as  y  increases. 

In  fact,  the  current  time  is  already  90%  of  t^  when 


(119) 


One  can  easily  show  that  if  the  assumption  a  /a  >>  1  is  removed,  the 

g  ° 

crack  growth  at  t  =  0  90  t^  will  be  even  less  than  given  in  Eq.  (119). 
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These  observations  are  important  because  they  tell  us  if  the  stress 
intensity  factor,  Eq.  (116),  is  valid  over  most  of  the  time  period 
0  <  t  <  t^  or  if  the  effect  of  finite  sheet  width  must  be 
included  when  predicting  failure  time.  For  three  typical  values  of 
m,  Eq.  (119)  yields 


a 

a 


10  ,  m  =  1 

3.3  ,  m  =  0.5 

[  1.6  ,  m  =  0.2 


(120) 


Surprisingly,  most  of  the  time  required  for  failure  is  consumed 
while  the  crack  is  still  relatively  small,  especially  for 
0  <  m  <  0.5. 

As  a  point  of  comparison,  if  the  bracketed  term  in  the  integrand 

in  Eq.  (117a)  is  replaced  by  (a  /a  Y)^m  for  all  0  <  t  <  tf,  the 

go  r 

failure  time  is  derived  by  setting  a/aQ  =  00  to  find 


2m  a*- 

t  =  _ B. 

f  (2+1/m) 
ir 


8T 


X  C_a 
m  2  o 


1/m 


-2 (1+1/m) 


(121) 


which  brings  out  very  simply  the  effect  of  stress  on  failure  time. 

This  is  the  same  result  as  would  be  obtained  by  neglecting 
2  2 

(N  /N  )  in  Eq.  (97)  at  the  outset  of  the  analysis, 
o  og  J 

Equations  (118)  and  (121)  will  be  compared  to  experimental  data 


on  Solithane  50/50,  which  is  a  crosslinked,  amorphous,  polyurethane 
rubber.  The  experimental  data,  which  are  shown  in  Fig.  9,  have  been 
normalized  with  respect  to  the  stress  oB,  where 
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/sr 

°®  ~  V  3ir 


(122) 


is  the  critical  stress  for  the  onset  of  crack  growth  in  an  elastic 

plate  having  Poisson's  ratio  v  *  1/2  and  Young's  modulus 
2 

E  ■  4(1  -  v  )/C  (°°)  ■  3/C  (®) ;  thus,  no  crack  growth  occurs  for 

oo  V  V 

2 

log(o/a<M)  <  0.  Also,  all  data  have  been  reduced  to  a  common 

temperature  of  0°C  by  recognizing  that  the  material  is  essentially 

thermorheologically  simple  [5];  according  to  Eq.  (103)  the  effect 

of  temperature  can  be  introduced  by  replacing  t^  with  t^/a,j,  I*1 

Eq.  (118).  Since  log(t/aT)  ■  log  t  -  log  aT»  the  theory,  if  valid, 

implies  experimental  failure  data  obtained  at  different  temperatures 

can  be  superposed  by  means  of  horizontal  translations  of  magnitude 

log  aT  so  as  to  form  a  single  curve;  the  failure  data  at  the 

temperatures  indicated  in  Fig.  9  were  shifted  to  the  data  at  0°C. 

(One  can  interpret  the  reduced  failure  data  as  being  the  results  of 

tests  conducted  entirely  at  0°C.)  We  should  add  that  the  creep 

compliance  in  Fig.  4  was  formed  in  the  same  way;  in  fact,  the  values 

of  log  a^,  used  to  shift  the  failure  data  are  those  obtained  by 

shifting  the  creep  compliance.  That  a  single  creep  compliance  curve 

was  obtained  is  a  check  on  the  thermorheological  simplicity 

assumption;  the  fact  that  the  failure  data  superpose  extremely  well 

(as  seen  in  Fig.  9)  using  the  same  values  of  log  a^  helps  to 

substantiate  the  assumption  that  T  and  o  are  constant,  and  checks  the 

m 

underlying  fracture  theory  itself. 
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In  view  of  the  power  law  fit  shown  in  Fig.  4,  we  set  m  *  0.5 
in  Eq.  (118)  and  carry  out  the  integration  to  find 


2  2 
0  o  I 

2  ml 


c  9 

( — 2_)2 

V2 


{! + 1 02  - 2  f +  ln 

v  O  c*  r\J 


(123) 


where  a^  depends  on  stress  through  Eq.  (117b).  When  ag/aQ  >  30  this 

equation  reduces  to  Eq.  (121)  after  setting  m  =  1/2.  The  latter 

equation  shows  the  theory  is  a  straight  line  with  slope  of  -1/3  when 
2 

log(o/o  )  is  plotted  against  log  t  ;  it  is  observed  in  Fig.  9 

that  this  is  indeed  the  slope  of  the  experimental  data  in  the  low 

stress  level  range  ((o/o  )^  <  10^'^). 

Referring  to  Eq.  (121).  we  see  that  the  influence  of  the  fracture 

properties  on  the  relation  between  a  and  tf  in  the  straight  line 

region  is  through  the  combination  (o  1,1).  Had  the  fracture  data 

m  1 

been  given  in  terms  of  a  in  [5],  rather  than  o/ob)  this  combination 

of  properties  could  have  been  found  by  matching  the  theory  and 

experiment  at  any  point  in  the  straight  line  range.  (Note  that 

o  I  and  T  appear  by  themselves  in  the  high  stress  level  range  of 

Eq.  (123);  if  the  theory  is  applicable  in  tnis  range  one  can  then 

determine  o  1^  and  T  separately.)  Knauss  [5]  obtained  the  fracture 

energy  itself  from  other  experiments  on  swollen  rubber,  in  which 

state  the  material  is  essentially  elastic,  and  then  used  the 

corresponding  value  of  to  normalize  the  failure  data.  He  found 
-2 

T  =  2.41  x  10  lb/in.,  which  corresponds  tc  =  8.40  psi. 

_2 

(Knauss  actually  reported  the  value  of  I'  =  3.21  x  10  lb/ir.  [5], 
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but  this  value  was  obtained  by  using  the  plane  stress  equation  for 

critical  stress  in  an  elastic  sheet;  later  it  will  be  found  that 

a  is  much  less  than  sheet  thickness,  which  means  the  plane  strain 

version  should  have  been  used.  In  order  to  correct  the  reported 

2 

T  value  we  have  multiplied  it  by  (1  -  v  )  *  3/4.  Now,  express 
Eq.  (121)  in  terms  of  (o/o  )  and  solve  for  o  1,  by  noting  that  the 
straight  line  in  Fig.  9  intercepts  the  stress  axis  at 
t^  =10®  sec  =  10® min;  hence 

o  I1  =  172  x  103  psi  .  (124a) 

By  taking  into  account  Eq.  (24b),  a  lower  bound  to  the  maximum 
failure  stress  is  obtained: 

om  >  86  x  103  psi  (124b) 

where  the  equality  is  used  if  is  constant  throughout  the  failure 

zone.  The  ideal  theoretical  strength  based  on  failure  of  a  regular 

arrangement  of  carbon  atoms  is  approximately  thirty-times  this 

value  [19];  however,  considering  the  irregular  nature  of  real 

networks,  the  presence  of  shear  forces  at  the  tip  (see  Fig.  3)  and 

that  Eq.  (124b)  is  the  lower  bound  and  not  necessarily  a  ,  this 

m 

numerical  result  does  not  seem  unreasonable. 

Inasmuch  as  both  o  I,  and  T  are  now  known,  we  can  plot  the 

ml  r 

remainder  of  the  solution  as  given  by  Eq .  (123).  While  having  the 

correct  shape,  the  theory  underpredicts  failure  time  in  the  high 

2 

stress  level  range  above  log  (a/am)  ~  0.8,  say;  this  stress 
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corresponds  to  an  overall  strain  of  approximately  5%.  Since 
tf  'v  nj*  for  m  *  1/2,  even  a  small  amount  of  nonlinearity,  such 
as  blunting  of  the  crack  tip  and  consequent  reduction  of  stress 
intensity  factor,  could  have  a  significant  effect.  In  fact,  we 
find  the  theory  agrees  with  all  of  the  experimental  results  if 
the  stress  intensity  factor  in  Eq.  (116)  is  replaced  with 

Nq  =  (1  +  2e)-1  a  (125) 


(It  is  to  be  noted  that  the  effective  modulus  of  the  sheet  is  a 
constant  430  psi  since  the  time  scale  for  crack  growth  in  Fig.  9 
is  far  beyond  the  viscoelastic  range  in  Fig.  4.) 

Even  though  the  stress  at  the  tip  is  very  high  (Eq.  (124b)), 
the  strain  is  not  because  practically  equal  triaxial  tension  exists. 
Using  a  typical  value  of  bulk  modulus,  2  x  10^  psi,  the  Young's 
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(ii) 


logfa/o^)2  =  1.6  : 


—8 

10  x  10  cm  ,  t 
32  x  10  ®  cm  ,  t 


0+ 

(126b) 

0.9  t 


where  Eq.  (125),  rather  than  Eq.  (116),  has  been  used  to  calculate 
stress  intensity  factor.  Although  tne  accuracy  of  the  theory  in 
predicting  failure  zone  lengths  comparable  to  interatomic  spacing 
(Eq.  (126a))  is  suspect,  it  is  believed  the  scale  is  at  least 
correct  in  view  of  the  agreement  between  experiment  and  theory  in 
the  low  stress  level  range;  if  a  were  constant,  for  example,  a  slope 
of  (-1/2),  rather  than  (-1/3),  would  be  predicted  for  the  graph  in 
Fig.  9.  (In  an  attempt  to  bracket  the  actual  failure  zone  size, 
stationary  and  moving  cracks  were  observed  with  the  TAMU  scanning 
electron  microscope.  Figure  11  shows  a  tip  which  was  moving  at 

-3 

approximately  10  in. /min  at  77°F.  Cracks  at  lOx  the  magnification 

of  this  figure  were  also  studied  visually  but  the  failure  zone  was 
still  too  small  to  be  observed;  because  of  excessive  discharging, 
with  and  without  a  coating,  the  resolution  of  the  crack  boundary 
was  not  as  good  as  in  Fig.  11.) 

The  instantaneous  value  of  a/a  is  found  from  Eqs.  (92b)  and 
(97)  to  be  a  decreasing  function  of  N  .  Therefore,  the  largest  a/a 
for  the  entire  set  of  experiments  occurs  at  t  =  0+  under  the  lowest 
stress.  We  find 


~  (max)  =  14  x  10  2  min  .  (127) 

a 

Referring  to  Fig.  4,  we  see  (a/ 3a) (max)  is  practically  within  the 


73 


range  of  the  generalized  power  law  creep  compliance;  since 
-4 

(a/a)  'v  Nq  and  Eq.  (127)  is  based  on  aQ  (rather  than  a),  use  of 
Eq.  (96)  with  m  =  1/2  is  justified  in  all  predictions  of  failure  times. 
Moreover,  the  maximum  fracture  initiation  time  is  approximately 
(a/3a)(max)  and  therefore  t^  is  negligible  relative  to  the  range  of 
tf  in  Fig.  9. 

Additionally,  as  required  by  the  underlying  theory,  the  change 

in  stress  intenity  factor  over  each  time  increment  a/a  is  small 

for  the  entire  experimental  range  of  behavior;  this  observation 

follows  from  Eq.  (127)  and  the  fact  that  the  increase  in  N  during 

o 

most  of  the  propagation  time  is  small  according  to  Eq.  (120). 


A  Long  Strip  under  Constant  Strain 


The  previously  determined  creep  compliance  and  fracture 
properties  of  Solitnane  50/50  will  be  used  in  conjunction  with 
Eq.  (97)  to  predict  constant  velocity  crack  propagation  in  the  long 
strip  shown  in  Fig.  10.  Provided  a  >  1.5b,  the  stress  intensity 
factor  in  a  strip  clamped  to  rigid  grips  is  [6], 


b(l-  vV1/2 

2tt 


(128) 


where  a  is  the  stress  which  must  be  applied  to  the  clamps  to  produce 

the  strain  e.  Although  this  applied  stress  and  strain  could  vary 

with  time,  the  experimental  results  are  for  the  case  in  which  the 

strain  is  held  constant  and  the  strip  is  in  a  fully  relaxed  state 

2 

outside  the  neighborhood  of  the  tip;  thus  a  =  E  e/(l  -  v  ),  which 
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when  substituted  into  Eq.  0  28)  yields 

No  *  E-(ii>1/2c  <129) 

where  we  have  set  v  =  1/2. 

Theoretical  upper  and  lower  bounds  on  strain  for  which  stable 
growth  exists  in  the  Solithane  are  obtained  by  using  Ec.  (94); 
there  results 


-1.89  <  log  e  <  -0.81 


(130) 


These  bounds  are  drawn  with  dashed  lines  in  Fig.  10. 

Except  at  high  strains,  the  predicted  velocities  are  seen  to 

be  in  very  good  agreement  with  the  experimental  data.  (If  the  stress 

intensity  factor  in  Eq.  (129)  is  divided  by  (1  +  2e),  as  in  Eq.  (125), 

the  predictions  at  the  high  strains  are  likewise  found  to  agree.) 

2  2 

Whenever  Nq/N  can  be  neglected  in  the  denominator  in  Eq.  (97) 

•  6 

we  find  a  *v  e  ,  and  therefore  the  theory  plots  as  straight  lines  in 
Fig.  10  with  slope  of  1/6;  the  upper  strain  limit  for  this  behavior 
is  found  to  be  at  log  e  =  -1.2. 

It  should  be  pointed  out  that  Mueller  and  Knauss  [6]  gave  a 
value  of  T  which  is  three  times  that  reported  in  Knauss'  constant 
stress  study  [5]  and  used  above.  This  larger  value  cannot  be 
reconciled  with  the  existence  of  failure  data  under  low  values  of 
conftant  stress  in  [5],  and  therefore  was  not  used  in  any  of  the 
calculations  made  here. 


For  the  range  of  experimental  data  in  Fig.  10,  Eq.  (92b)  with 
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corrected  Tor  finite  strain  as  mentioned  above,  yields 

0.70  x  10“8  <  a  <  53  x  10"8  cm  .  (131)  j 

The  largest  value  of  a/3a,  after  converting  to  T  *  0°C,  is 

_3 

0.9  x  10  min;  this  result  implies,  of  course,  that  the  generalized 
power  law  in  Fig.  4  is  valid  for  the  entire  experimental  range. 

Variable  Loading  and  Properties;  the  Linear  Cumulative  Damage 

Rule 

Some  consequences  of  Eq .  (95)  concerning  crack  growth  and  failure 
will  be  examined  and  then  compared  with  experimental  results  on  two 
different  particulate  composites:  solid  propellant  and  asphaltic 
concrete.  For  reasons  of  clarity,  we  first  introduce  some  simplifying 
assumptions.  Although  they  may  not  be  satisfied  exactly  in  practice, 
it  is  believed  there  are  many  applications  in  which  the  resulting 
error  will  have  little  effect  on  the  total  time  required  for  failure. 

Except  for  some  later  observations,  we  assume  (i)  n  is  constant. 
The  other  material  properties  (C^,  T,  o^I^)  may  be  functions  of 
time  (which,  for  example,  could  be  due  to  a  transient  temperature 
and/or  aging)  but  we  assume  (ii)  their  variation  during  most 
intervals  a/a  is  small.  Also,  motivated  by  Eq.  (115b),  we  assume 
(iii)  fracture  initiation  time  for  the  initial  loading  (as  well  as 
for  subsequent  reloading  if  there  are  periods  for  which  a  =  0)  is 
negligible.  Finally,  it  is  assumed  ( iv)  that  the  stress  intensity 


factor  depends  on  only  one  length  parameter  defining  either  the 
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crack  size  "a"  (say)  or  local  geometry  of  the  body  "h"  (say). 

From  dimensional  considerations, 

N  =  /a  L  (132) 

o 

or 

N  =  v'h  L  (133) 

o 

where  "L"  is  a  linear  function  of  the  stresses  applied  to  the  body; 
e.g.,  see  Eqs.  (116)  and  (128).  The  size  "a"  may  represent  the 
radius  of  an  isolated  internal  (penny-shaped)  crack  or  the  half-length 
of  an  isolated  through-crack.  The  parameter  "h"  could  be  the 
separation  of  two  rigid  boundaries,  where  the  crack  length  or 
diameter  has  to  be  at  least  somewhat  greater  than  h;  recall  that 
a  >  1.5b  is  sufficient  for.  the  plane  problem  in  Fig.  10. 

It  should  be  emphasized  that  if  th  "*e  are  periods  when  the 
stress  intensity  factor  is  negative  (due  to,  for  example, 
cyclic  loading  having  a  zero  mean  value)  the  analysis  is  not 
invalidated.  For  the  above  assumptions  (ii)  and  (iii)  together  with 
Eq.  (95)  imply  a  during  growth  periods  depends  on  the  instantaneous 
stress  intensity  factor  and  not  on  either  its  history  or  the  stress 
history.  The  following  equations  will  be  given  without  explicitly 
writing  a  =  0  whenever  L  <  0;  however,  one  should  set  L  =  0  during 
all  periods  in  which  L  is  actually  negative. 

Supposing  first  that  Eq.  (132)  applies,  integration  of  Eq.  (95) 
and  setting  a  =  00  yields  an  implicit  equation  for  failure  time,  t^, 
in  terms  of  the  initial  half-length,  a^,  and  material  property  and 
load  histories, 
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1  =  - 


,  (2n+l)l 

a  a  it'  / 

1/n 

o  n 

4 

*  , 

Cx  1/n  2(l+l/n) 

[— ]  - 2 — 2 —  dt 


(134) 


a_  i; 

m  1 


This  equation  can  be  expressea  in  a  more  recognizable  form  by 


introducing  a  new  quantity,  t£,  where 


t  =  t  (t)  = 
c  c 


2  2 
no  I, 
m  1 


l2(1+1/r) 


4r 


.  (2n+l) 

a  \  it  C, 

on  1 


1/n 


(135) 


Then  Eq.  (134)  becomes  simply 


(136) 


It  is  easily  demonstrated  by  means  of  Eq.  (134)  that  t^(t)  would 
be  the  failure  time  for  the  body  if  it  had  timewise  constant  applied 
loads  and  constant  material  properties  equal  to  those  that  exist 
in  the  actual  variable  loading  problem  at  the  current  time  t. 

If  L  and  the  properties  are  piecewise  constant  in  time, 
then  Eq.  (136)  becomes 
N  At. 

1  “  Zf  T1  (137) 

j=l  cj 

where  At^  is  one  of  the  time  periods  over  which  the  body  has  constant 
applied  loads  and  properties,  t  would  be  the  failure  time  for  the 
body  if  it  had  these  timewise  constant  loads  and  properties  from 
t  =  0  to  failure,  and  is  the  number  of  steps  that  are  needed  to 
produce  failure  in  the  actual  variable  loading  case. 

Equations  (136)  and  (137)  are  two  forms  of  the  well-known 
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linear  cumulative  damage  rule  which  has  been  used  by  engineers  to 
predict  failure  times  of  specimens  and  structures  (usually  with 
acceptable  accuracy)  subjected  to  varying  loads,  given  experimental 
data  on  the  failure  time  t  for  constant  loads.  This  rule,  but  in 
a  slightly  different  form,  was  originally  proposed  by  Palmgren  [20] 
and  used  by  Miner  [21]  as  an  empirical  method  of  predicting  failure 
under  variable  amplitude  cyclic  loading  of  metals;  it  is  often  called 
"Miner's  Hypothesis"  or  "Miner's  Law." 

Another  feature  of  Eq.  (134)  may  be  brought  out  by  predicting 

failure  under  constant  loading  rate,  L.  For  simplicity  we  assume 

the  mechanical  and  fracture  properties  are  constant.  Substituting 

L  =  L  t  in  Eq.  (134),  and  defining  the  "critical  stress  intensity 

factor"  N  by 
cr  1 

N  =  /T~  L  t  (138) 

cr  ox  x 


(which  is  the  stress  intensity  factor  existing  at  the  time  of  failure, 

tx>  say,  referred  to  the  initial  crack  size  a  ,  say)  yields 
r 


N  =  i 
cr 


,,  ^  3/2  2  2  .* 

(3n  +  2)  a  o  I  L 
ox  m  1 


4r 


i  J2n+1>r 

V  ci 


1/n 


n/(3n+2) 

.  (139) 


The  insensitivity  of  to  loading  rate,  L,  when  0  <  n  <  1  is 
especially  to  be  noted. 

Now,  return  to  Eq.  (95),  which  is  not  limited  to  a  stress 

intensity  factor  of  the  form  in  Eq.  (132).  Still  assuming  the 

fracture  properties  are  constant  we  m3y  express  these  properties 

in  terms  of  t  ,  N  ,  and  a  ;  Eq.  (95)  becomes 
x  cr  ox 
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da 

dt 


,3n  +  2v  .3qX'  ,No  ,2(l+l/n) 

'  aT  Mt  HN  ' 

T  x  cr 


(140) 


This  equation  enables  one  to  predict  crack  growth  and  failure  of  a 
structure,  given  the  critical  stress  intensity  factor  determined 
experimentally  in  a  simple  specimen  having  initial  crack  length 
2a^x  and  failing  in  time  t^.  In  writing  Eq.  (140)  we  have  assumed 
thermorheological  simplicity  and  selected  a^,  to  be  unity  at  the 
temperature  of  the  experiment  used  to  establish  Ncr*  The  value  of 
=  a^,(T)  shown  in  Eq.  (140)  corresponds  to  the  instantaneous 
temperature  of  the  structure  of  interest. 

If  instead  of  taking  a  specimen  to  failure,  one  simply  measures 
the  crack  velocity,  a x,  corresponding  to  an  instantaneous  stress 
intensity  factor,  Nx  (which  need  not  be  restricted  to  the  forms  in 
Eqs.  (132)  and  (133)) ,  then  the  crack  growth  equation  (95)  becomes 

da  =  (V(l+l/n)  (141) 

An  advantage  of  this  experimental  method  over  that  leading  to  Eq.  (140) 
is  nonlinear  effects  can  be  minimized  in  the  test  by  measuring 
velocity  early  in  the  growth  period.  However,  Eq.  (140)  does  not 
require  a  velocity  measurement  and,  therefore,  may  be  more  desirable 
if  nonlinearity  Js  not  a  problem. 

Consider  now  sttess  intensity  factor  Eq.  (133),  in  which  h  is 
constant.  Let  us  define  failure  in  this  case  as  the  time  at  which 
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the  crack  length  reaches  some  preselected  value,  a^,  say.  Then 
all  of  the  results  Eqs.  (134)  -  (137)  apply  directly  to  this 


case  after  making  the  substitution 
1/n 


o 

n 


1+1 /n 


a  -  a 
mx  o 


(142) 


in  Eqs.  (134)  and  (135). 

The  sensitivity  of  crack  velocity,  and  therefore  failure  time, 

to  stress  intensity  factor  is  clearly  revealed  by  Eqs.  (140)  and 

(141).  Suppose,  for  example,  Eq.  (140)  is  used  to  predict  failure 

time  of  a  structure  having  n  =  0.2,  a  constant  value  of  Nq,  initial 

flaw  length  a  *  a  ,  and  failure  is  defined  to  occur  when  the  crack 
o  ox 

grows  to  ten  times  its  initial  length.  Suppose  further  that  the 

experimental  failure  time  is  t  *  1  minute  and  for  the  structure 

x 

a_  *  1  and  N  *  N  /2.  We  find  the  structure  fails  at  tr  =  10  days. 
T  o  cr  f  1 

If,  however,  Nq  =  N^/4,  t^ien  ~  100  years!  Obviously,  the 
ability  to  predict  service  life  of  viscoelastic  structures  is  very 
much  dependent  on  the  accuracy  of  the  stress  intensity  factors. 
However,  because  of  the  variability  of  material  properties  and 
environmental  and  loading  factors,  one  can  expect  considerable 
scatter  in  failure  times  when  a  large  number  of  structures  are 
involved,  which  is  already  a  well-known  empirical  fact. 

Before  turning  to  the  application  of  the  theory  to  specific 
materials,  it  ahould  be  emphasized  that  the  cumulative  damage  rule 
is  not  generally  valid  in  the  range  for  which  n  is  time-dependent. 
A  simple  illustration  is  provided  by  Eq.  (97)  using  stress  intensity 
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factor  Eq.  (132).  Even  with  a  constant  value  for  m,  the  presence 
2 

of  Nq  in  the  denominator  precludes  writing  the  right-hand  side  in 
the  form  f(a)  x  g(L),  which  is  necessar/  if  the  linear  cumulative 
damage  rule  is  to  be  derived. 

Failure  of  Solid  Propellant;  Failure  data  on  several  different 
solid  propellants  (which  consist  of  a  cross-linked  rubber  matrix 
filled  with  approximately  75  vol  %  hard  particles  having  diameters 
typically  on  the  order  of  10-200  microns)  under  constant  uniaxial 
stress  have  been  reported  by  Bills  et  al.  [22,23]  to  fit  the  equation 


where 


iiil  =  (J^k 
Do  V 


(143) 


oT  *  o(l  +  e)  =  true  stress 
o  =  engineering  stress  (constant  during  test) 

£  =  engineering  strain  at  failure 
o  -  engineering  failure  stress  below  glass  transition 

o 

temperature 
*  =  tc/aTq 

g 

q  =  constant  =  10 

a^,  =  temperature-dependent  time-scale  factor 

t  =  failure  time 
c 

D (ip )  =  uniaxial  propellant  creep  compliance  (=  D  ij>n  in 

K 

the  range  of  <p  for  which  the  experimental  failure 


data  exist,  where  D  and  n  are  essentially 

K 

independent  of  temperature) 
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and 

D  -  D(o)  . 
o 

Also,  1.5  <  k  <  2.5  is  approximately  constant  for  any  one  propellant. 
In  the  referenced  work  [22],  the  relaxation  modulus,  E(i|>) ,  is  used 
instead  of  D(if));  to  obtain  Eq.  (143)  we  have  applied  the  quasi-elastic 
approximation 

E(*)  D«»)  <*  1  (144) 


which  is  quite  good  for  the  usual  range  of  n,  0.15  <  n  <  0.25  [24]. 
It  was  found  that  failure  did  not  occur  when  the  tensile  stress  was 
below  some  relatively  small  value,  which  evidently  corresponds  to 
the  stress  for  which  Nq  <  N  (see  Eq.  (94b)).  In  writing  Eq.  (143) 
it  is  implied  that  the  stress  is  above  this  lower  limit. 

The  theoretical  equation  (135)  will  be  cast  in  the  same  form  as 
Eq.  (143);  in  doing  this  we  obtain  a  partial  check  on  the  theory  and 
enable  the  empirical  constants  (k,q)  to  be  evaluated  in  terms  of 
basic  material  properties.  First,  note  that  Eq,  (48),  together  with 
the  fact  that  v  -  1/2  and  that  propellants  are  (approximately) 
thermorheologically  simple, 

D 

C  (O  =  3  D(0  =  3  Dn5n  =  3  —  tn  (145) 

v  k  n 

aT 

and  therefore 


(146) 
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Next,  introduce  the  stress  intensity  factor  for  an  (internal) 
penny-shaped  crack  of  radius  a  [25] 

,,  -  Z2a 

Nop  =  —  o  (147) 

and  define  the  ratio  R  as 

R  =  actual  stress  intensity  factor  of  the  most  rapidly 
growing  crack  in  specimen/N^ 

Then  substitute  L  =  R  N^/Za  into  Eq.  (135),  assume  cr,_,  -  a  and  find 
that  it  becomes  identical  to  Eq.  (143)  if 

k  =  2(1  +  n)  (148a) 


and 


n  2  l  ,  «• 

’  '  2  'l  lcgRI 


3a  A  D 
o  n  o 


(148b) 


In  order  to  interpret  the  parameters  in  Eq.  (135)  or  (148)  in 
terms  of  the  properties  of  solid  propellant  and  its  constituents,  it 
is  essential  to  consider  the  size  of  the  failure  zone  relative  to 


the  particle  sizes  and  spacing.  Figure  12  shows  this  zone  at  the 
tip  of  a  through-crack  in  a  sheet  of  propellant.  It  is  seen  that 
the  zone  is  large  relative  to  the  particles,  and  consequently  the 
overall  or  "effective"  properties  of  the  composite  propellant 
are  the  appropriate  ones  to  use  in  the  fracture  model;  on  the  other 
hand,  if  the  failure  zone  had  been  entirely  within  a  relatively 
large  volume  of  matrix  material,  the  matrix  mechanical  and  failure 


«*>»*•  *«•«*{»  wu- 
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properties  would  be  the  ones  to  use  in  the  fracture  theory.  Now, 
empirical  Eq.  (143)  was  established  using  specimens  wichout  cracks 
except  for  naturally  existing  flaws;  the  fact  that  the  correlation 
of  stress  and  failure  time  for  all  propellants  studied  is  in  terms 
of  propellant  compliance  (rather  than  that  of  the  matrix)  is  believed 
to  imply  that  the  crack  tip  geometry  is  similar  to  that  seen  in 
Fig.  12,  at  least  during  most  of  the  time  required  for  failure. 

As  further  evidence  of  this  point,  one  can  usually  see  very  fine 
surface  cracks  (with  sizes  on  the  order  of  the  largest  filler 
particles)  in  a  specimen  under  tensile  stress;  and  these  cracks  grow 
only  a  small  amount  during  most  of  the  time  required  for  failure  [24]. 
Since  propellant  contains  a  large  percentage  of  much  smaller  particles, 
it  is  tentatively  concluded  that  typically  2a^  >  200  microns,  and 
that  the  failure  zone  is  significantly  larger  than  most  of  the 
particles;  since  the  compliance  of  the  surrounding  continuum  is 
essentially  that  of  the  propellant,  n  =  0.2,  and  consistent  with 
observations  we  predict  from  Eq.  (120)  that  gross  failure  follows 
shortly  after  the  largest  crack  has  grown  approximately  60%  in  size. 

It  is  to  be  observed  that  Eq.  (148a)  predicts  k  =  2.4  when 
n  =  0.2,  which  is  on  the  high  side  of  the  experimentally  determined 
range  1.5  <  k  <  2.5.  That  the  reported  values  of  k  are  not  larger 
is  believed  to  be  primarily  due  to  the  r-iuction  of  the  effective 
stress  intensity  factor  with  strain,  as  suggested  previously  in  the 
study  of  Soli  thane. 


As  further  evidence  of  the  validity  of  the  theory, 
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Bills  et  al.  [22,23]  have  applied  the  linear  cumulative  damage  rule, 
Eq.  (136),  together  with  Eq.  (143),  to  the  failure  prediction  of 
solid  propellant  specimens  and  structures  under  monotonic  and 
cyclic  temperatures  and  stresses  with  considerable  success. 

One  final  observation  concerns  a  method  of  predicting  fracture 

initiation  in  propellant  structures  that  was  originally  suggested 

and  successfully  applied  by  Swanson  [26];  it  was  motivated  by  a 

technique  used  widely  for  metals.  Essentially  the  method  consists 

of  applying  a  constant  strain  rate  to  a  precracked  specimen  and 

noting  the  critical  stress  intensity  factor  and  time  when  the 

specimen  breaks.  It  is  then  assumed  that  a  crack  in  a  more  complex 

structure  will  start  to  propagate  when  the  stress  intensity  factor 

reaches  the  value  determined  from  the  laboratory  specimens,  even 

though  the  stress  history  in  the  structure  may  be  different.  By 

referring  to  Eq.  (140)  or  (141)  we  can  see  why  such  a  method  is 

successful,  at  least  for  some  stress  histories.  Inasmuch  as 

n  'v  0.2,  the  exponent  2(1  +  1/n)  'v  12,  which  implies  there  is  very 

slow  crack  growth  (relative  to  laboratory  specimen  time  scale  and 

crack  size)  when  <  N  ,  while  the  velocity  is  relatively  high 

when  N  >  N  .  However,  one  must  be  aware  of  the  limitations  of 
o  cr 

this  method  as,  for  example,  structural  failure  can  occur  even 

though  the  stress  intensity  factor  is  always  less  than  N  ;  of 

course,  failure  time  might  be  very  long  relative  to  the  experimental 

time  for  failure,  t  . 

x 

Fatigue  Crack  Growth  in  Asi’.altic  Concrete;  As  a  result  of  a 


large  amount  of  fatigue  testing  it  has  been  found  that  the  crack 
growth  per  cycle,  da/dN,  can  be  expressed  as  the  simple  power  law 
[27-29], 


da 

dN 


B  N 


om 


(149) 


where  B  and  3  are  constant  under  fixed  environmental  conditions 
and  input  wave  shape,  and  N^m  is  the  amplitude  of  an  oscillating 
stress  intensity  factor.  For  beams  and  plates  on  foundations  with 
sufficiently  high  rigidity  8  =  4  at  77°F  [27,28].  Probably  due  to 
strain  effects  (e.g.,  crack  tip  blunting),  values  of  8  somewhat 
less  than  4  have  been  reported  for  unsupported  specimens  [27,29]. 
Also,  it  has  been  noted  that  8  tends  to  increase  with  decreasing 
temperature  [29]. 

We  shall  show  that  the  theoretical  fracture  model,  Eq.  (95), 
is  consistent  with  these  findings.  First,  introduce  the  function 
w  =  w(t)  which  defines  the  wave  shape  of  the  stress  intensity  factor. 


N 

o 


om 


where  the  maximum  value  of  stress  intensity  factor  during  a  cycle, 

which  is  N  ,  may  vary  from  cycle-to-cycle .  We  make  the  reasonable 
om 

assumption  that  crack  growth  per  cycle  is  small,  which  implies  w(t) 
is  essentially  equal  to  the  wave  shape  of  the  externally  appli*'! 
loading.  Now  substitute  definition  Eq.  (150)  into  Eq.  (95)  and 
integrate  over  a  cycle  having  period  t^;  writing  the  crack  growth 
per  cycle,  Aa,  as  da/dN  we  find 


1 


da  _  2(l+l/n) 

dN  -  Bt  Nom 


(151) 


where 


rc,X  ^(20+1)1 1/n r t+tp 

1  n _  w  _ _ 

4  1/n  2  2 

it  F  °m  X1 


(152) 


Note  that  F  and  ate  shown  in  the  integrand  as  they  could  be 

time-dependent,  depending  on  the  complexity  of  the  wave  shape  and 
whether  or  not  the  crack  is  stationary  during  part  of  the  cycle; 
the  exponent  n  and  are  assumed  constant  during  each  cycle. 

Equality  of  empirical  Eq.  (149)  and  theory  Eq .  (151)  requires 


B  =  B  and  6  =  2(1  +  1/n) 


(153) 


Interpretation  of  these  results  turns  out  to  require  that  we 
differentiate  between  a  failure  zone  which  is  entirely  within  the 
bitumen  matrix  ("micro-zone")  and  one  which  encompasses  many  aggregate 
particles  ("macro-zone") ,  with  the  latter  situation  being  analogous 
to  that  in  Fig.  12  for  solid  propellant.  For  a  macro-zone  the 
properties  in  Eq,  (152)  and  the  exponent  (n)  are  those  of  the 
asphaltic  concrete  mposite;  while  there  may  be  exceptions,  the 
typical  range  of  the  exponent  for  the  composite  is  0  <  n  <  0.5, 
at  xeast  for  temperatures  not  exceeding  100°F,  with  n  decreasing 
monotonically  with  decreasing  temperature  until  it  is  approximately 
zero  at  temperatures  below  the  glass  transition  temperature  [e.g.,  30,31] 


On  the  basis  of  the  latter  observation,  we  conclude  6  >  6  for  a 
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macro-zone. 

This  value  of  8  is  considerably  larger  than  the  values 
(6  4)  commonly  reported  for  asphaltic  concrete  above  40°F.  Thus, 

we  must  consider  the  possibility  of  the  failure  zone  being  entirely 
within  the  bitumen.  Pursuing  this  question  we  found  from  microscope 
studies  at  room  temperature  that  the  failure  zone  is  not  at  all 
like  that  in  Fig,  12.  Rather,  the  bitumen  matrix  is  so  weak  (it  is  a 
non-crosslinked  polymer  having  a  much  smaller  molecular  weight  than 
the  propellant  binder)  that  the  failure  zone  of  a  single  crack 
appears  to  be  confined  between  individual  particles  instead  of 
being  spread  out  over  many.  If  there  is  a  sufficiently  large 
neighborhood  of  bitumen  around  the  tip  (which  is  certainly  possible 
if  a  is  comparable  to  the  zone  size  in  solithane,  for  example) 
then  the  value  of  n  and  other  material  properties  in  Bt  , 

Eq.  (152),  would  he  those  for  the  bitumen.  Except  at  temperatures 
close  to  and  below  the  glass-transition  temperature  of  the  bitumen, 
the  value  n  =  1,  which  implies  from  Eq.  (153)  that  8  =  4,  is  believed 
to  be  appropriate.  Besides  predicting  the  experimentally  observed 
value  of  8  at  small  strains,  it  is  to  be  noted  that  the  creep 

compliance  of  the  bitumen  is  g even  by  D(t)  =  except  at  short 

times  or  at  low  temperatures  [32];  in  addition,  a  state  of  high 

triaxial  stress, which  exists  at  the  tip  according  to  the  earlier 

analysis  and  in  the  neighborhood  of  a  tip  confined  between  rigid 
particles,  is  known  to  suppress  the  glass  transition  temperature  of 
polymers  [33]  and  thereby  extend  the  temperature  range  for  which 
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8  =  4  is  predicted. 

Carrying  this  interpretation  somewhat  further,  let  us  consider 
the  micro-zone  to  be  a  precursor  of  a  macroscopic  crack  having  total 
length  much  larger  than  a  and  particle  spacing.  Then  since  the 
local  amplitude  of  stress  intensity  factor,  N  ,  is  proportional  to 
the  far-field  stress  amplitude  normal  to  the  precursor,  o^,  say,  we 
may  write 


N 

om 


0-54) 


1/2 

In  the  simplest  situation  Eq.  (128)  implies  g  =  (3b/8n)  '  , 

P 

where  2b  is  the  distance  between  adjacent  particles  and  we  have 
set  v  =  1/2.  More  generally,  g^  can  be  expected  to  depend  on  particle 
shape  and  size  as  well  as  spacing.  Because  of  the  relative  rigidity 
of  the  aggregate  particles  compared  to  the  bitumen,  dimensional 
considerations  imply  g^  is  essentially  independent  of  modulus;  then 
assuming  v  =  1/2  for  the  bitumen, the  correspondence  principle  [13] 
implies  Eq .  (154)  is  valid  for  viscoelastic  media.  Now,  for  a 
sufficiently  long  crack. 


o'  =  N’  /A,  (155) 

y  om  1 

where  N'  is  the  stress  intensity  factor  one  would  calculate  for  the 
om 

asphaltic  concrete  when  viewed  as  a  homogeneous  continuum,  and 
at  least  somewhat  larger  than  the  particle  spacing  in  the 
neighborhood  of  the  crack  tip.  Substituting  Eq .  (155)  into  (154) 
and  the  result  into  Eq .  (151),  we  obtain  3  result  which  is  formally 
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identical  to  original  Eq.  (151);  viz., 


—  =  B* (N*  )2(1+1/n) 
dN  V  om; 


(156) 


where  depends  on  particle  spacing  and  geometry,  in  addition  tc 

bitumen  properties,  while  N'  is  the  "effective"  stress  intensity 

om 

factor  for  the  composite  material.  In  general,  both  B’  and  N' 

t  ora 

could  be  functions  of  the  number  of  cycles  N. 

While  the  above  description  of  crack  growth  in  asphaltic 
concrete  represents  at  best  a  tentative  model  and  a  partial  check 
on  the  theory,  it  does  bring  out  the  importance  of  the  microstructure 
(i.e.,  particles)  in  the  fracture  of  composite  materials.  Proper 
interpretation  of  experimental  results  obviously  requires  that  the 
failure  zone  size  be  known  relative  to  microstructural  dimensions, 
and  that  the  possibility  of  its  size  changing  appreciably  with 
temperature  be  recognized. 

With  n  constant,  and  N'  having  the  form  of  either  Eq.  (132) 

om 

or  (133),  cumulative  damage  rules  analogous  to  Eqs.  (136)  and  (137) 
are  easily  derived.  The  only  difference  in  the  result  is  that  t  is 
replaced  by  the  number  of  cycles,  N,  and  t  by  ,  where  is 
the  number  of  cycles  required  to  fail  a  specimen  under  constant 
amplitude,  N^,  and  constant  B^.  The  rule  Eq.  (137)  then  becomes 
identical  to  Miner’s  Hypothesis. 
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5.  Related  Theories 

As  mentioned  in  the  Introduction,  the  existing  theories  of 
crack  propagation  in  viscoelastic  media  have  already  been  reviewed 
in  Knauss'  recent  article  [2]  and,  therefore,  such  an  effort  will  not 
be  repeated  here.  Rather,  we  shall  simply  compare  some  results  of 
this  paper  with  closely  related  ones  established  by  others. 

The  governing  equation  derived  by  Knauss  [5]  fcr  growth  of  a 
central  crack  in  a  large  sheet  (see  Fig.  9)  and  one  derived  by 
Mueller  and  Knauss  [6]  for  a  cracked  strip  (see  Fig.  10)  have  the 
same  form  as  Eq.  (91).  Namely,  the  relations  are  essentially  those 
for  an  elastic  material  except  creep  compliance  replaces  the  elastic 
compliance.  Although  there  is  some  difference  in  certain  details 
it  is  surprising  how  similar  the  results  are,  considering  the  fact 
that  the  method  of  derivation  is  so  different  in  each  case. 

In  [5],  the  work  done  by  the  continuum  on  the  failing  material  is 
approximated  using  the  same  relation  as  for  an  elastic  material 
(viz.,  1/2  force  x  displacement)  along  with  some  other  analytical 
approximations  at  the  crack  tip.  The  result  is 

D(a/a)  =  ~~  .  (157) 

7io  a 

In  contrast,  the  approach  in  [6]  to  predicting  crack  propagation  in 
a  long  strip  is  to  let  the  crack  grow  in  steps.  With  the  crack 
stationary,  the  stress  in  the  failure  zone  is  assumed  to  decay 
linearly  with  time  until  it  vanishes;  the  tip  is  then  instantaneously 
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advanced  a  distance  a  and  the  stress  again  decays  linearly  in  time. 
With  v  *  1/2  the  governing  equation  is 

G(a/a) - ~  (158) 

2be  E 

00 

where  G  is  the  difference  of  two  integrals  involving  the  uniaxial 
creep  compliance;  for  Solithane  it  is  shown  in  [6]  that 
G(a/a)  =  D(n/3a).  Both  Eqs.  (157)  and  (158)  were  derived  using  a 
singular  stress  distribution  along  with  the  assumptions  of  plane 
stress  and  a  constant  value  of  a.  As  we  have  already  discussed, 
when  a  is  much  less  than  sheet  thickness  (it  is  for  Solithane) 
plane  strain  theory  should  be  used,  which  accounts  for  one  difference 
between  the  results  in  [5,6]  and  the  theory  in  this  paper. 
Furthermore,  a  is  not  really  constant;  referring  to  Figs.  9  and  10, 
for  example,  with  constant  a  the  slopes  of  the  straight  portion  of 
the  curves  are  predicted  to  be  (-1/2)  and  (1/4),  as  compared  to  the 
actual  values  of  (-1/3)  and  (1/6),  respectively.  Finally,  since 
x|/n  s  1/3  for  all  0  <  n  <  1,  according  to  Fig.  6,  the  theory  in 
this  paper,  Eqs.  (91)  and  (92a),  shows  that  the  use  of  D(a/3a) 
on  the  left-hand  side  of  Eq.  (158)  is  essentially  correct,  while 
the  argument  a/a  in  Eq.  (157)  is  three-times  too  large.  In  later 
work  dealing  with  rapidly  changing  stresses,  Knauss  and  Dietmann  [7] 
removed  the  assumption  of  plane  stress  and  considered  relatively 
general  crack  geometries,  but  the  other  features  in  [6]  (i.e., 
stepwise  crack  propagation,  singular  stress,  and  linear  decay  with 
time  of  stresses  near  the  tip)  were  retained. 
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Turning  to  the  prediction  of  fracture  initiation  time,  t^, 
Williams  [3]  used  a  spherical  flaw  geometry  in  order  to  establish 
this  time.  For  a  constant  stress  applied  at  t  =  0,  he  obtained  a 
relation  which  is  very  close  to  an  equation  for  critical  stress  in 
an  elastic  body.  Translating  Williams'  solution  to  the  notation  in 

Eq.  (108)  for  plane  strain,  we  find  his  result  is  obtained  if 

/ 

C^(t^)  is  replaced  by  2Cv(t^)  -  C^(o).  Thus,  the  two  theories  agree 
qualitatively  for  t^  >  0  and  predict  the  same  critical  stress 
intensity  factor  when  t  =  0.  Predictions  for  time-varying  stresses 
do  not  lend  themselves  to  such  a  direct  comparison.  Additionally, 
when  the  stress  intensity  factor  is  specialized  to  that  for  a 
penny-shaped  crack,  Eq.  (108)  predicts  the  same  fracture  initiation 
time  as  derived  by  Wnuk  and  Knauss  [4]  for  a  viscoelastic-perfectly 
plastic  solid  having  a  constant  yield  stress. 

6.  Concluding  Remarks 

A  simple  theory  of  crack  propagation  has  been  developed 
and  successfully  applied  to  three  different  unfilled  and  filled 
polymeric  ma.'_erials  under  constant  and  variable  loading.  Although 
the  continuum  was  assumed  linearly  viscoelastic,  the  nature  of  the 
failure  zone  is  quite  arbitrary  and,  therefore,  could  include 
material  which  is  highly  nonlinear,  rate  dependent,  and  even 
discontinuous.  Under  mildly  restrictive  conditions,  an  equation  for 
fracture  initiation  time  was  derived  and  found  to  be  very  similar 
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to  the  elasticity  relation  for  critical  stress,  except  a  "secant 
compliance"  appears  in  place  of  an  elastic  constant. 

Some  influence  of  finite  strains  on  crack  velocity  was  noted 
for  the  three  materials  studied.  The  general  tendency  appears 
to  be  that  the  actual  crack  velocity  falls  below  the  linear  theory 
prediction  as  the  strain  is  increased;  it  was  suggested  that  an 
important  cause  is  geometric  nonlinearity,  in  that  crack-tip  blunting 
at  high  strains  could  reduce  the  effective  stress  intensity  factor. 

An  ad  hoc  modification  for  finite  strain  effects  in  Solithane,  in 
which  a  simple  product  form  was  used,  brought  the  theoretical  and 
experimental  results  together  at  high  strains.  As  to  whether  or  not 
one  can  normally  expect  to  be  able  to  account  for  geometric  and 
material  nonlinearities  in  the  continuum  through  such  a  simple 
modification  is  not  known  at  this  time.  A  point  of  encouragement  is 
that  the  relaxation  modulus  of  many  polymers  can  be  expressed  as 
a  product  function  of  strain  and  time  [3], 

Although  only  the  opening  mode  has  been  analyzed  in  this  paper, 

it  seems  reasonable  to  expect  the  governing  equations  for  fracture 

initiation  time  and  crack  growth  in  shearing  modes  to  be  analogous 

to  those  developed  here  and,  therefore,  to  be  similar  to  equations 

for  critical  stress  intensity  factors  in  elastic  media.  Thus,  for 

the  skew-symmetric  mode  [34]  (i  ^  0  along  plane  of  crack  prolonga- 

xy 

tion)  it  is  conjectured  that,  in  analogy  with  propagation  Eq.  (91), 

t 


Cv<‘a>  =  ^2 
tiN4 
2 


(159) 
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where  energy  T  is  not  necessarily  equal  to  that  for  the  opening  mode, 

and  the  stress  intensity  factor  N2  is  defined  in  the  singular 

shear  stress  distribution  near  the  tip,  x  =  N„//cT;  while  for 

r  xy  2  1 

the  antiplane  mode  [34]  (x  i  0  along  the  plane  of  crack  prolongation), 

zy 

2J(t  )  «  ~  (160) 

where  J(t)  is  the  creep  compliance  in  shear  and,  just  as  for  K^, 
is  defined  in  the  singular  distribution  t  =  N^/v^".  Equations 
for  predicting  a  in  these  two  modes  probably  are  similar  to  Eq.  (92b), 
but  determination  of  their  exact  form  requires  further  analysis. 

When  two  or  three  stress  intensity  factors  exist  simultaneously, 
and  it  is  assumed  that  the  fracture  energy  is  the  same  at.  in  the 
opening  mode,  we  add  the  work  input  to  the  failure  zone  for  each 
mode  and  find  that  Eq.  (91)  is  to  be  replaced  by 

(N2  +  N2)  C  (t  )  +  2N2  J(t  )  -  —  (161) 

o  2  v  a  3  a  it 

whenever  Nq  >  0.  The  opening  mode  does  not  contribute  energy  if 
Nq  <  0  since  the  displacement  in  the  y-direction  is  then  zero; 
we  must  therefore  set  Nq  ;  0  in  Eq.  (161)  unless  Nq  is  positive. 

Without  further  analysis  it  is  not  clear  what  the  length  of 
the  failure  zone  should  be  under  combined  loading.  Nevertheless,  by 
means  of  the  following  idealized  problem  the  nature  of  the  interaction 
can  be  demonstrated:  Assume  N_=  0,  N  >0,  the  failure  stress 
=  om  ■  constant  with  the  x-  stress  component  negligible,  and  a 
maximum  principal  stress  criterion  in  which  0^  is  the  same  regardless 
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of  the  relative  magnitudes  of  Nq  and  Nj.  By  analogy  with  the  deriva¬ 
tion  of  Eq.  (92b)  we  deduce 


a 


N 


o,2 


ifcr-r  +  N;r~> 


(162) 


The  major  problem  in  generalizing  the  opening  mode  theory  to 
include  shearing  modes  is  believed  to  be  that  of  predicting  the 
direction  of  crack  tip  motion,  as  it  is  very  likely  that  the  crack 
will  not  extend  in  its  original  plane.  An  example  of  the  complexity 
of  the  problem  may  be  seen  in  Knauss*  paper  [35],  in  which  a  crack 
was  forced  to  propagate  in  Solithane  in  the  antiplane  shear  mode. 
Crack  growth  was  found  to  occur  by  the  opening  of  semi-penny-shaped 
cracks  which  straddle  the  main  crack  front  at  an  angle  of  45  degrees; 
i.e.,  the  individual  planes  of  crack  extension  have  normals  parallel 
to  the  direction  of  maximum  principal  stress.  In  the  more  general 
case  when  two  or  three  modes  of  loading  exist  together,  one  can 
expect  the  direction  of  crack  motion  to  depend  at  least  on  some 
average  of  the  maximum  principal  stress  directions  over  the  failure 
zone.  That  the  propagation  direction  does  not  depend  solely  on  the 
state  of  stress  at  the  crack  tip  (£^  =  0)  may  be  argued  by 
considering  the  stresses  in  rubber,  in  which  case  v  =  1/2;  since 
the  opening  mode  stresses  at  the  tip  are  essentially  equal,  the 
direction  of  propagation  when  and/or  are  not  zero  would  be 
relatively  insensitive  to  Nq,  which  contradicts  observed  behavior. 

Finally,  let  us  briefly  consider  the  problem  of  opening-mode 
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adhesive  fracture.  The  critical  stress  required  to  cause  propagation 
of  an  isolated  debond  between  a  linear  elastic,  incompressible 
(v  =  1/2)  continuum  and  a  rigid  substrate  is  identical  to  that  for 
cohesive  fracture  of  the  continuum  itself  except  the  adhesive 
fracture  energy  r  say,  replaces  the  cohesive  fracture  energy  [36]. 
Moreover,  it  is  a  simple  matter  to  show  that  the  extended 
correspondence  principle,  which  was  used  in  the  derivation  of 
Eq.  (91),  is  valid  as  long  as  the  substrate  is  rigid  and  the 
continuum  is  homogeneous  and  incompressible.  If  these  restrictions 
are  essentially  met,  we  suggest  that  Eq.  '1)  can  be  used  in  predicting 
propagation  of  the  debond  crack  after  T  is  replaced  by  T  . 

au 

However,  it  is  anticipated  that  some  modification  of  Eq.  (92b)  for 
failure  zone  length  will  be  required  as  both  shear  and  normal 
stresses  act  along  the  interface. 
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We  first  examine  the  order  of  the  integral  AI,  Eq.  (31),  when  the 


failure  stress-difference  is  given  by 


A afU)  *  A  Cq 


(A-l) 


in  the  interval  0  <  C  <  a^,  where  0  <  <  a,  A  and  q  are  non-zero 
constants,  and  0  <  q  <  1/2;  for  this  range  of  q  Eq.  (15)  is  met  but 
not  Eq.  (17).  The  integral  in  Eq.  (31)  is  rewritten  as  follows: 


(A-2) 


3/2 

The  second  integral  is  readily  shown  to  be  of  order  £  by  expanding 
the  logarithm  in  a  power  series  in  the  variable  /£/£'  and  letting 
£  -*■  0.  The  first  integral  is  expressed  in  terms  of  a  new  variable  of 


■a 

ri 

= 

+ 

o 

integration,  y  =  /£ ’/£  ,  and  together  with  Eq.  (A-l)  it  becomes 

} 


f  1 

r/ot  /C 

-  2A  «1+<! 

o 

y21  1 

2  -  y  In 

Y  +  1 

Y  l 

0 

y  -  1 

dy 


(A-3) 


Again  divide  the  range  of  integration  into  two  parts, 


ra  1 


2A  K 


1+q 


{  *  } 


(A-4) 
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Since  we  are  interested  in  the  order  of  Eq.  (A-4)  as  £  -*  0,  the  upper 
limit  and  the  constant  are  taken  to  satisfy  the  inequality 
0  «  y^  <  this  choice  of  together  with  the  series  expansion 


for  y  »  1, 


y  In  1  -  I  =  2  +  -2-  +  _S_  +  . 
r  •  1  3y2  5/ 


renders  Eq.  (A-4)  as 


(A-5) 


=  2A  51+q  |b  + 


3(1  -  2q) 


5(3  -  2q) 


12  '  q 


(A-6) 


where  q  ^  and  B  is  non-zero  and  is  independent  of  £.  When 
1 

q"2’ 


-  2A  51*1  {b  +  i  in  i-  +  •••} 

*  n 


(A-7) 


Equations  (A-6)  and  (A-7)  together  with  Eq.  (A-2)  imply  that 


AI  =  0(C1+<1) 


(A-8) 


when  0  <  q  <  —  and 


AI  =  Oa3/2  In  O 


(A-9) 


when  q  =  y.  It  therefore  follows  from  Eqs.  (27)  and  (30)  that  the 
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order  of  v  is  also  given  by  Eqs.  (A-8)  and  (A-9)  for  0  <  q  <  y  and 
q  *  respectively.  Note  that  if  q  -*■  0,  the  displacement  v  0(5) 
which  implies  the  crack  cross-section  approaches  a  straight-sided 
wedge  very  close  to  the  tip. 

The  limiting  case  q  =  0  corresponds  to  a  constant  stress  over 
0  <  5  <  a^j  in  order  that  a ^  be  continuous  to  the  left  at  5  “  0,  in 

3/2 

accordance  with  Eq.  (15),  we  must  take  A  =  0.  For  this  case  v  =  0(5  )• 

Prediction  of  the  behavior  of  the  stresses  near  the  tip  is 
accomplished  in  a  similar  manner.  First,  we  record  the  resultant 
stresses,  Eq.  (19b),  after  substituting  Eq.  (28)  for  o^: 


o  =  a  =  o  (0)+cr,  (0)  [—  tan'1,/!-  -  1] 
x  y  f  f  V  5, 


Aaf(5) 
^(5  +  5L) 


(A-10) 


where  5^  >  0. 

1/2 

The  second  term  is  of  order  5-^  ,  which  follows  from  a  power 

series  expansion  of  the  arctangent  in  the  variable  /5^/a.  Given 
the  power  law,  Eq.  (A-l)  ,  the  order  of  the  integral  in  the  third 
term  is  established  by  following  essentially  the  same  step3  used  to 
analyze  AI  (i.e.,  Eqs.  (A-2)  -  (A-7)).  We  find 


1  ^  (O 

-  sc  — 5 —  dc 

*  ‘U^v 


7  Ba  EJ  +  °«l/2>  =0^4 


7  cj/2  In  ^  +  0(^/2> 


(A-ll) 


where  is  non-zero  and  is  independent  of  E,^.  Thus, 


=  o  =  of (0)  +  0(£?) 


o  =o 
x  y 


(A-12) 


when  0  <  q  <  —  and 


Ox  =  Oy  =  of (0)  +  0(cj'  In  Cl) 


(A-13) 


when  q  =  y.  Note  that  Eqs.  (A-12),  (A-13),  and  (19a)  predict  these 
stresses  to  have  a  vertical  tangent  at  =  0;  it  can  be  shown  that 
only  for  the  unlikely  situation  in  which  of(£)  has  a  horizontal 
tangent  at  £  =  0  and  the  integral  1^,  Eq.  (39),  vanishes  will  the 
stresses  °x(^)  anc^  ay(4^)  also  have  a  horizontal  tangent  at  £^. 

The  above  order  properties  for  stress  and  displacement  are 
different  from  those  in  [37}  because  do^/d£  was  assumed 
to  be  bounded  in  [37], 


B.  Criterion  for  Crack  Growth  when  o  •+  00 

m 


At  the  outset  the  normalized  failure  stress  distribution, 
f  =  0f/°m»  *-s  assumed  to  satisfy  only  the  following  physically-based 


assumptions : 


(i)  f(£)  is  piecewise  continuous  and  non-negative 


(ii)  lim  f(£)  =  f(0)  >  0 
£  0+ 
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Now,  let  us  show  that  a  vanishes.  To  this  end,  write  Eq.  (22) 
in  the  fora 


it  N 

\/a  I - 2. 

1  o„ 


(B-l) 


m 


Supposing  that  the  stress  intensity  factor  is  finite,  we  obtain 


Sa  I,  0  as  o 
1  m 


Moreover,  the  integral  1^,  Eq.  (23),  cannot  vanish  on  account  of 

the  above  two  conditions  on  f(C);  hence,  a  -+-  0  as  a  Letting 

m 

a  -»■  0  in  Eq.  (23)  yields  1^  =  2  f (0)  and,  in  turn,  Eq.  (B-l)  yields 


lira  [2VoT  a  f  (0) / tt  ]  =  N 
v  m  o 

a  °° 

m 


(B-2) 


With  continuous  crack  growth,  da/dt  >  0,  the  total  time  needed 

for  the  crack  to  propagate  the  distance  a,  viz.,  -  t^,  approaches 

zero.  Therefore,  assuming  f(0),  o  ,  and  N  are  continuous  functions 

ra  o 

or  time,  or  else  are  independent  of  time,  Eq.  (B-2)  shows  that  a  is 
essentially  independent  of  time  during  the  period  t^  <  t  < 

In  view  of  these  results,  we  will  evaluate  the  displacement, 

Eq.  (41),  and  then  the  fracture  energy,  Eq.  (57),  for  the  case  of  a 
very  small,  constant  a.  First,  write  displacement  Eq.  (41)  by 
making  an  integration-by-parts  and  the  following  change  of  variables: 


£  V  f 

n  =  -  ,  n  =  —  ,  f  =  — 

a  la  a 

m 


(B— 3) 


where  0  <  (n , n ,f )  <  1 


The  result  is 


106 


2n  v 


0  a  v 
m 


C  (0) 


f  (otn,)[2(-5-)1/2  -  In 
1  nl 


+  */n" 


-  /n 


dn. 


ft' 


cv(c  “  T) 


'  o 


f  (an,) [2(^-)1/2 
1  nl 


-  In 


+  ^n". 


-  /r) 
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(B-4) 


In  which  <  t  < 


n  =  [a(t)  -  x]/ot 


(B-5) 


in  the  first  line,  and 


n  -  [a0v)  -  x]/o  (B-6) 

in  the  second  and  third  lines.  Also 

C'v( t)  -  d  cv/dt  (B-7) 

is  finite  and  positive  for  all  t  >  0  for  real  materials;  it  is  a 
simple  matter  to  show  that  Cv(t)  is  the  creep  compliance  of  a  plate 
of  infinite  extent  in  the  z-direction  (i.e.,  -  0)  and  subjected 

to  a  constant  stress  in  the  x-  or  y-direction. 

Equation  (B-4)  yields 


11m  [2v  vAy,]  *  Cv(0)  f  (0)  [2^  -  (1  -  n)  ln  L±A i 
*  “  1  -  ^ 


(B-8) 


The  second  and  third  lines  in  Eq.  (B-4)  vanish  on  account  of  the 
finiteness  of  C’  and  /*,  and  the  fact  that  tn^sa^  co;  reCall 
that  t1  <  t  <  t2  where  (t2  -  tj)  f  0  as  »  for  da/dt  >  0. 

Now,  change  the  variables  in  fracture  energy  Eq.  (57)  using 
definitions  in  Eq .  (B-3)  for  r,  and  f: 


F  =  °m  f(cxr>>  f*  • 


(B-9) 


Rewrite  this  equation  to  read 


tlm)  h  * 

m 


(B-10) 


and  let  affl  ->  °°  (which  implies  a  ->  0+)  : 


lim  t271  r/om  "  f(0)  ~  ■  lim  dn  . 

a  ■>  on  a  a J I 
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o  m 


(B-ll) 


The  process  of  taking  the  limit  under  the  integral  and  inside  the 
derivative  is  valid  on  account  of  the  piecewise  continuity  of  f(£) 
for  0  <  4  <  a  and  the  boundedness  and  continuity  of  the  derivative 
in  Eq.  (B-10)  with  respect  to  0  <  n  <  1  and  a  +  0+.  Integrate 
Eq.  (B-ll)  and  then  use  Eq.  (B-8)  to  find 
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lim  [ 2n  r/a2  a]  -  2 C  fO)  f2(0)  .  (B-12) 

m  v 

a  ->•  °° 
m 

The  quantity  /a  o^  f(0)  can  be  eliminated  between  Eqs.  (B-2)  and 
(B-12)  to  obtain  the  equation  for  critical  stress  intensity  factor 
of  a  brittle  elastic  solid, 
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(B-13) 


The  initial  value  theorem  for  Laplace  transforms  [17]  may  be  used 
in  conjunction  with  Eq.  (45)  to  find 


cv(0) 


Ml  -  v2(0)i 

E(0) 


(B-14) 


which  when  substituted  into  Eq.  (B-13)  yields  Eq.  (59). 

Graham  [25]  used  a  global  energy  criterion  and  the  singular 
stress  distribution  for  a  large  centrally-cracked  sheet  and  for  a 
large  body  with  a  penny-shaped  crack.  He  obtained  results  for  the 
critical  stress  which  are  identical  to  those  obtained  from  Eq.  (B-13) 
when  Nq  is  specialized  to  these  two  geometries. 


C.  Work  Input  to  Failure  Zone  for  a  One-Term  Representation 
of  v 


Work  input  based  on  approximate  displacement  Eq.  (72a)  will 
be  derived  and  compared  with  energy  Eq.  (91).  The  following  two 
classes  of  failure  stress  are  assumed: 
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(C-l) 


(ii) 


Of  =  < 


m 


m  v, 
b 


;  0  <  v  <  v, 

-  -  b 


v  >  v. 


(C-2) 


where  the  quantities  o-t  V  v  a,  and  b,  as  well  as  ^  and  n  in 
power  law  Eq.  (66),  are  assumed  to  be  positive  constants;  in  reality 
they  may  depend  on  tip  velocity,  but  for  the  problem  at  hand  it  is 
not  necessary  to  show  this  dependence  explicitly. 

The  second  distribution  rigorously  implies  a  =  »,  and  we  shall 
use  this  limit  in  evaluating  the  integrals  A  However,  from  a 
Practical  standpoint,  if  b  is  sufficiently  large  the  influence  of 
the  stress  of  outside  of  some  finite  distance  a^,  say,  will  be 
negligible;  it  is  assumed  that  «ef/0  «  1,  where  6  is  the  length 
parameter  appearing  in  Eqs.  (2)  -  (4). 

The  approximate  work  input  derived  using  displacement 
Eq.  (72a)  and  the  above  failure  stresses  will  be  denoted  by  T 

ap 

We  ratio  the  results  to  f  as  given  by  Eq.  (91), 


r  =  ?  V  V  No 


(C-3) 


and  obtain,  finally: 
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where  b  >  1.  These  ratios  are  plotted  in  Fig.  C-l  for  a  few  different 
values  of  the  constants.  The  ratio  for  the  case  a  *  1/2  is  rather 
large  at  n  =  0  because  this  exponent  is  close  to  the  value 
(a  *  1/3)  for  which  the  integral  Eq.  (39),  diverges.  We  see  that 
for  most  of  the  predictions  the  error  is  no  greater  than  33%. 
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Figure  4.  Creep  Compliance  in  Uniaxial  Tension  for  Solithane  50/50 
(experimental  data, - >  after  Mueller  and  Knausa  [6]) 
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Figure  7.  Normalized  Displacement  Used  in  Energy 
Integrals  T  and  r. 
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Figure  11.  Crack  in  Solithane  50/50. 


